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ABSTRACT 


This  paper  investigates  the  effect  of  cooling  a  hot 
Planckian  emitter  upon  its  fluence  spectrum  .  A  sequence  of 
models  of  increasing  complexity  is  developed  to  determine 
the  effects  of  various  aspects  of  cooling  upon  the  spectrum, 
such  as  initial  temperature,  density,  and  ionization  state 
of  the  plasma.  Spectra  are  calculated  for  radiating  plasmas 
composed  of  different  atomic  number  materials  (carbon,  alu¬ 
minum,  copper,  and  plutonium)  at  initial  temperatures  of 
0.02  -  10  KeV,  and  initial  densities  of  lE-t-25  -  1E+29 
atoms/m3,  to  observe  the  effects  of  these  parameters  on  the 
fluence  spectrum.  The  change  in  material  and  binding  energy 
for  some  spectra  at  the  low  energy  end  produces  a  second, 
prominent  but  smaller  peak.  The  resulting  non-Planck ian 
spectra  can  be  approximated  with  two  or  more  Planckian  basis 
functions  having  different  temperatures. 
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TIME- ! NTLGRATED  SPECTRUM 
OF  A  RADIATIVELY  COOLING  PLANCK  I  AN  EMI  TTt.fi 

I .  Introduct ion 

This  paper  deals  with  computing  the  spectral  distribu¬ 
tion  of  thermal  radiation  emanating  from  an  extremely  hot 
body,  and  how  cooling  effects  its  spectrum  taken  over  time. 

Background 

Astronomers  have  long  been  interested  in  the  spectral 
character  of  radiation  coming  from  hot  bodies  such  as  stars, 
quasars,  and  supernovas,  to  obtain  clues  to  their  makeup, 
temperature,  age,  and  internal  processes  that  occur  in  their 
life  cycle.  It  is  generally  recognized  that  nuclear  fusion 
is  the  process  that  heats  the  cores  of  stars  to  extremely 
high  temperatures  (millions  of  degrees  Kelvin),  and  causes 
them  to  radiate  enormous  quantities  of  energy. 

Currently,  there  is  much  research  effort  going  on  world¬ 
wide  to  achieve  a  controlled  and  sustained  source  of  fusion 
energy  as  an  alternative  to  other  sources  of  energy.  Fusion 
experiments,  then,  can  serve  as  likely  sources  for  studying 
the  spectral  nature  of  hot  bodies  (plasmas)  undergoing  spon¬ 
taneous  cooling,  once  the  power  source  is  shut  off. 
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One  example  of  an  experimental  setup  that  is  capable  of 
achieving  high  temperatures  in  a  plasma  (although  not  yet  to 
the  point  of  fusion  ignition)  is  the  imploding  liner  pinch 
illustrated  in  Figure  1,  adapted  from  (1:315). 


Fig.  1.  Imploding  Liner  Pinch 


A  plasma  is  contained  within  a  cylindrical  thin  metal  tube 
or  liner  through  which  a  large  electric  current,  Ij  ,  is 
passed.  The  ax i a  1 -d i r ec ted  current  induces  an  azimuthal 
magnetic  field,  B«  ,  around  the  liner.  Both  current  and 
magnetic  field  interact  together  to  produce  a  radial, 
inward-directed  force,  F.  ,  that  causes  the  liner  and 
plasma  to  implode.  The  compressive  force  and  inertial 
momentum  of  the  liner  are  eventually  halted  by  the  increas¬ 
ing  plasma  pressure  from  within,  but  by  that  time  compres- 
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sion  has  superheated  both  plasma  and  liner  past  the  point  of 
vaporizing  and  ionizing  the  liner  material  which  becomes  a 
part  of  the  plasma. 

In  order  to  predict  the  spectral  behavior  of  a  hot  body, 
such  as  the  plasma  in  the  imploding  liner  pinch,  through  a 
mathematical  model,  it  is  sometimes  assumed  that  the  temper¬ 
ature  of  material  and  radiation  is  the  same,  therefore  the 
distribution  of  radiant  energy  emitted  from  the  body's 
surface  can  be  described  as  Planckian.  This  study  also 
assumes  this.  However,  integrating  the  spectrum  over  time, 
i.e.,  the  time-integrated  spectrum  or  fluence  spectrum, 
should  show  the  behavior  to  be  non-Planck ian. 

Prob lem 

The  purpose  of  this  study  is  to  investigate  the  effects 
of  cooling  on  the  fluence  spectrum  of  a  Planckian  emitter. 

In  addition,  the  spectrum  is  compared  to  a  Planckian  curve- 
fit  spectrum  to  determine  the  validity  of  approximating  the 
fluence  spectrum  with  Planckian  basis  functions. 

Scope 

The  study  is  computat ional  and  analytic  rather  than 
experimental,  so  the  results  are  predictive  in  nature. 

Three  parameters  are  varied  as  initial  conditions  for  the 
Planckian  emitter  to  observe  their  effects  during  cooling. 
These  are  temperature,  density,  and  material  composition. 

The  initial  temperature  range  is  0.02  -  10  KeV,  and  the 
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emitter  is  allowed  to  cool  to  about  1/50, 000th  of  the  ini¬ 
tial  temperature.  Densities  range  from  1E+25  -  1E+29  atoms 
per  cubic  meter,  and  material  composition  is  one  of  the 
following  elements:  carbon,  aluminum,  copper,  or  plutonium. 

Approach  and  Assumptions 

A  sequence  of  models  is  considered  that  increases  in 
complexity  according  to  the  type  of  internal  energy  assumed 
to  be  contained  within  the  emitter,  whether  it  be  material 
energy,  radiation  energy,  or  a  combination  of  both,  etc. 
These  models  and  their  assumptions  are  as  follows: 

Model  1  -  Radiative  Cooling  with  Constant.  Material  Heat 
Capacity,  but  No  Radiation  Heat  Capacity. 

-  The  volume  is  constant;  this  model  excludes  cooling  by 
expansion. 

-  Local  thermodynamic  equilibrium  exists  between  matter 
and  radiation. 

-  The  temperature  is  uniform  throughout  the  volume  at 
any  instant,  but  varies  with  time. 

-  The  heat  capacity  remains  constant  throughout  the 
cooling  process. 

Model  E  -  Radiative  Cooling  with  Temperature  Consistent. 
Radiation  Heat  Capacity,  but  No  Material  Heat  Capacity. 

-  The  above  assumptions  apply  except  the  radiation  heat 
capacity  varies  with  temperature. 
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Model  3  -  Radiative  Cooling  with  Constant.  Material  Heat 
Capacity  and  Temperature  Consistent.  Radiation  Heat 
Capac i tv . 

-  The  assumptions  of  Model  1  and  S  taken  together  apply 
to  this  model,  with  clarification  that  the  heat  capacity 
contributed  by  material  is  constant,  and  the  heat  capac¬ 
ity  contributed  by  radiation  changes  with  temperature. 
Model  4  -  Radiative  Cooling  with  Temperature  Consistent 
Heat  Capacity. 

-  Model  3  assumptions  are  retained  except: 

-  The  material  heat  capacity  is  allowed  to  vary  accord¬ 
ing  to  the  ionization  state  of  the  atoms  which  is  depen¬ 
dent  on  temperature. 

Presentat ion 

Chapter  II  presents  the  general  analytic  form  and  theory 
behind  calculating  the  time-integrated  spectrum,  and  Chapter 
III  shows  the  numerical  approximation  for  the  analytic  form. 
Then  Chapters  IV  -  VII  present  the  models  in  order  as  dis¬ 
cussed  above  along  with  results  and  discussion  of  the  com- 
|kj  l.pd  spectra.  Conclusions  and  recommendations  are  included 
in  Chapter  VIII. 
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I  I .  Theory 


Let  us  start  with  a  definition  of  the  radiatively  cool¬ 
ing  Planckian  emitter  and  its  surroundings,  followed  by  some 
assumptions  about  how  they  interact  through  the  law  of 
energy  conservat i on ,  then  develop  the  concept  of  a  time- 
integrated  spectrum,  first  as  a  function  of  time,  then  as  a 
function  of  temperature.  From  the  latter  development,  it 
becomes  apparent  that  heat  capacity  plays  a  major  role  in 
the  spectrum  behavior. 

Definitions 

The  initially  hot  body  under  consideration  can  be  any 
collection  of  ions,  electrons,  and  neutral  atoms  from  a 
single  element  constrained  to  an  arbitrary,  constant,  finite 
volume.  For  initial  temperatures  of  0.02  to  10  KeV  under¬ 
taken  in  this  study  (1  KeV  =  11,600,000  K)  and  densities  of 
1E+25  to  1E+29  particles/m3,  an  optically  thick  plasma  is  an 
appropriate  description  for  the  hot  body.  Within  the  vol¬ 
ume,  the  material  particles  collide  with  one  another  giving 
rise  to  radiation  particles  or  photons  which  in  turn  are 
absorbed  or  scattered  by  the  material  particles.  At  any 
given  time,  a  fraction  of  these  photons  reach  the  volume 
boundary  and  radiate  outward,  away  from  the  body,  depleting 
the  plasma  of  its  internal  energy.  The  loss  of  energy  is 
equivalent  to  lowering  the  temperature  of  the  plasma,  hence 
it  cools  by  radiative  emission. 
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Besides  the  process  of  radiative  cooling,  two  other 
mechanisms  exist  for  cooling  the  body.  They  are  the  trans¬ 
fer  of  kinetic  energy  from  material  particles  to  the  wall  or 
force  field  defining  the  volume  boundary,  and  expansion  of 
the  volume;  however,  neither  cooling  mechanism  is  considered 
in  this  study. 

At  any  instant  in  time,  assuming  the  radiation  tempera¬ 
ture  is  the  same  as  the  kinetic  temperature  of  the  material 
particles  in  a  condition  of  local  thermodynamic  equilibrium, 
the  photon  energy  distribution  can  be  described  by  the 
Planck  function  at  some  average  temperature.  Thus,  the 
cooling  emitter  is  Planckian  in  nature. 

Just  beyond  the  plasma  boundary  lies  an  infinite,  empty 
universe  at  temperature  r*Q  Kelvin  into  which  the  radi¬ 
ation  escapes.  It  is  assumed  that  there  is  no  flux  of  radi¬ 
ation  into  the  volume. 

Energy  Balance 

Having  defined  the  cooling  emitter  and  its  surroundings, 
let  us  now  account  for  the  gain  and  loss  in  energy  of  the 
plasma  through  a  differential  energy  balance  equation.  In 
general,  the  change  in  total  internal  energy  of  the  plasma 
with  respect  to  time  is  equal  to  the  heat  flux  in  minus  the 
heat  flux  out  of  the  plasma,  plus  heat  coming  from  any 
internal  source  minus  heat  lost  to  any  internal  sink,  minus 
work  done  by  the  expanding  material.  This  is  expressed  in 
the  equation: 
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d  U  /  d  t  -  Q  flux  (a  ~  Q  finx  oat  +  Q  toarco  *3  wok  ~  ^  by  (  1  ) 

where  the  dot  over  a  variable  denotes  differentiation  with 
respect  to  time.  Since  there  is  no  flux  of  radiation  into 
the  body  from  the  outside  universe,  the  term  for  heat  flux 
in  is  zero.  Regarding  the  source  term,  the  supply  of  energy 
to  heat  the  plasma  has  stopped  (by  assumption)  so  this  term 
is  also  zero.  The  heat  sink  term  represents  energy  absorbed 
by  the  material  particles  in  endothermic  reactions,  but  it 
is  assumed  that  the  particles  do  not  react,  therefore  this 
term  vanishes.  One  more  term  vanishes  because  the  volume 
remains  constant  with  time:  the  work  term,  which  can  be 
defined  as  pressure  times  the  change  in  volume  with  respect 
to  time. 

Now  the  change  in  internal  energy  per  unit  time  is  that 
which  radiates  out  through  the  surface  area  in  the  same 
amount  of  time: 

di//d  t  =  t  =  -oT'.A  (2) 

where  crTj  is  the  energy  per  unit  area  emitted  per  unit  time 

from  the  total  surface  area  A  ;  T,  is  the  surface  tempera¬ 
ture;  and  0  is  the  Stef an-Bo 1 t zmann  constant.  In  the  above 
equation  the  minus  sign  signifies  a  loss  of  energy  to  the 
surrounding  universe.  The  term  <JT}  is  a  result  of  Stefan's 
law  for  black  body  radiation  (2:4),  and  is  equivalent  to  the 
Planck  distribution  of  energy  per  unit  area  per  unit  time 
Da(£,T,)  integrated  over  all  photon  energies  E*Kv  (  h.  is 
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Planck's  constant  and  v  is  frequency)  at  temperature  r,. 

This  surface  temperature  is  the  same  as  the  average  internal 
temperature  at  time  t,  and  is  hereby  designated  T(0« 


Concept  of  the  T i me- I ntearated  Spectrum 
Let 

f  r*  t  \  /  r*  t*  /  /  >  >  ^  r  «r  \ 

pac.c;  ■  U;;  -  >1  ,) 

be  the  Planck  function  that  describes  the  energy  distribu¬ 
tion  per  unit  area  per  unit  time  so  that 

[\A(E.t)AdE  =  <»[  T(t)]*A  (4) 

J  0 


If  the  hot  body  starts  with  initial  internal  energy  U 1  and 
is  allowed  to  cool  forever,  it  will  eventually  lose  all  its 
energy  such  that  the  following  expression  holds  true: 

U\  =  f  a[T(t)]'Adt  C5) 

J  o 

hence 


F 


(E.t)AdE 


d  t 


U  i 


(6) 


and  defining 


f(£)  -  f\A(E.t)Adt  (7) 

Jo 

equation  (7)  is  the  distribution  function  that  describes  the 
t ime- integrated  spectrum,  which  upon  normalizing  to  unit 
area  is  given  by 


S(£) 


f  /  r  > 


LkAi 


F 


C8) 
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Change  of  Variables  from  Time  to  Temperature 


So  far,  the  development  of  the  t i me- i ntegr ated  spectrum 
has  dealt  explicitly  with  the  time  behavior  of  the  radiation 
energy  distribution,  with  temperature  being  a  function  of 
time.  This  requires  knowledge  of  how  the  temperature  varies 
with  time,  which  may  not  be  well  known.  An  equivalent  and 
easier  approach  is  to  determine  the  behavior  as  a  function 
of  temperature  so  that  the  integration  of  equation  (7)  is 
performed  over  temperature  instead  of  time.  This  approach 
eliminates  the  need  to  know  the  time  behavior  explicitly, 
and  the  end  result  is  the  same:  the  fluence  spectrum  (or 
time-integrated  spectrum). 

Referring  back  to  equation  (7),  the  integration  variable 
is  changed  from  time  to  temperature: 

f(£)  -  fT,p/l(F.T)A~dT  (9) 

Jt  i  d  T 


where  the  initial  temperature  T ,  corresponds  to  t,  =» 0  ,  and 
the  final  temperature  T{  corresponds  to  rf»«®  which  may  be 
considered  as  T,= 0  ,  and 


Pi( E-T)A  = 


2n 


c  2  ft. 3 


(10) 


is  the  Planck  distribution  of  energy  per  unit  time.  Equa¬ 
tion  (10)  is  adapted  from  (2:24)  in  terms  of  energy  instead 
of  frequency.  A  is  Planck's  constant,  c  is  the  speed  of 
light  in  a  vacuum,  and  k  is  Boltzmann's  constant. 
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Before  proceeding  further,  dt/dT  in  equation  (9)  can  be 


transformed  by  considering  the  change  in  total  internal 
energy  as  a  total  derivative  with  respect  to  time,  volume, 
and  temperature  in  which  the  concept  of  heat  capacity  is 
introduced : 

d jJ  =  +  +  iVfdT) 

d  t  *  w{dtl  +  aTldtJ  C 


In  the  first  term  on  the  right  hand  side,  dU / dt  is  zero 
because  the  internal  energy  does  not  change  with  just  the 
passage  of  time  and  no  other  influences.  In  the  second 
term,  there  is  no  change  in  volume  V  as  stated  before,  so 
this  term  is  zero.  Only  the  third  term  remains  which  con¬ 
tains  within  it  the  definition  of  heat  capacity,  that  is 


Uf  i  \r 


/  1  />  \ 


;  •  neat  capacity  =  au/oi 
which  is  the  change  in  energy  per  unit  change  in  tempera¬ 
ture.  Combining  equations  (2),  (11),  and  (12)  and  rearrang¬ 
ing  terms  results  in  the  desired  transformation  for  : 

di/di  s  -M.il/flJ  /I  0  i3j 


Now  equation  (9)  becomes 


HE) 


CmnAjE.T)A 

aT'A 


d  T 


(14) 


Let  be  normalized  to  unit  area  in  the  following 

manner  with  T  as  a  parameter; 


1  1 


HElT) 


(15) 


P *(£»?)  / 
<jT4  i4 


2*  r  1 

!/i3(or4)Le£f/tr5-  1 J  uo; 


2nBfc4 

15c2/;3 


^(F;T)  = 


Jmp(E\T)dE  - 


15 

(nkT) 


.[-EL. 

4|  e(//*n_ 


(2:24)  (17) 


(18) 


(19) 


Incorporating  pCFjT")  into  equation  (14)  and  interchanging  the 


order  of  integration  yields: 


f(F)  -  fT,p(E;T)C(T)dT 

•  o 


(20) 


This  expression  says  that  the  time-integrated  spectrum  at 
photon  energy  F  has  a  value  equal  to  the  sum  of  the  released 
energy  with  a  certain  distribution  at  each  temperature  T 
from  the  initial  temperature  down  to  zero.  By  integrating 
equation  (20)  over  all  energy  F  and  dividing  the  result  into 
equation  (20),  the  result  is  the  normalized  fluence  spectrum 
equivalent  to  equation  (8)  at  any  energy  E  : 


fi(F) 


f  (F)  =  f (£) 

j*f(E)dE  F 


(21) 


Heat  Capacity 

Up  to  this  point,  we  have  explicitly  defined  the  normal¬ 
ized  Planck  function  {>( E',T)  in  equation  (18),  but  have  not 
explicitly  defined  the  heat  capacity  CC7")  needed  to 
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calculate  the  time-integrated  spectrum.  Since  heat  capacity 
is  the  change  in  internal  energy  with  respect  to  tempera¬ 
ture,  a  question  to  ask  is  how  the  internal  energy  of  the 
plasma  varies  with  temperature.  A  general  formula  for  the 
total  internal  energy  U  as  a  function  of  temperature  is: 

U(J)  =  UU(J)  *  Ut(J)  *  U m(T)  (22) 

where  d/m  is  the  material  energy;  d/f  is  the  internal  radi¬ 
ation  energy;  and  Uw  is  the  energy  due  to  latent  heats  of 
transformation.  is  comprised  of  several  terms  that 

describe  the  energy  required  to  change  the  material  state: 
that  from  solid  to  liquid  (latent  heat  of  fusion),  from  liq¬ 
uid  to  vapor  (latent  heat  of  vaporization),  and  once  in  the 
vapor  state,  to  ionize  or  strip  off  electrons  from  the  atoms 
(latent  heat  of  ionization).  Specific  forms  of  CCD  are 
deferred  to  the  corresponding  chapters  that  talk  about  how 
each  term  in  equation  (22)  affects  the  fluence  spectrum 
through  CCT")  . 
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Numer ical  Approx imat i on 


III  . 

Equations  (SO)  and  (SI)  ,  representing  the  non- 
normalized  and  normalized  fluence  spectrum,  respectively, 
are  difficult  to  solve  analytically.  Therefore,  numerical 
methods  are  employed  with  the  use  of  a  computer  to  approxi¬ 
mate  the  solutions  to  these  equations. 

For  convenience,  temperature  is  redefined  as  oskT  , 
having  units  of  energy.  Both  y  and  E  are  expressed  in  KeV 
throughout  the  remainder  of  the  paper. 

Simpson's  rule  of  integration  is  used  for  both  integra¬ 
tion  over  temperature  and  energy,  with  the  limits  of  inte¬ 
gration  altered  slightly  from  the  original  limits.  For  the 
lower  limits,  y  and  E  cannot  assume  the  value  of  zero 
because  either  condition  causes  division  by  zero  in  equa¬ 
tions  (20)  and  (21).  Therefore,  the  lowest  value  for  y  is  a 
final  temperature  of  yr>0  ,  and  that  for  E  is  E,„>Q 
An  upper  limit  of  E = »  is  likewise  impossible  to  evaluate 
numerically  in  a  computer  solution,  thus  in,**®  •  In  fact, 

£ki(k  does  not  need  to  be  more  than  about  3O0j  because  the 
upper  tail  of  the  spectrum  (plot  of  f( E)  or  $(E)  versus  S') 
always  decreases  monotonical ly  to  zero  by  E  =  *>,  for  all  forms 
of  C(0)  considered  in  this  study,  and  the  area  under  the  tail 
in  this  region  is  negligible  compared  to  the  area  beneath 
the  spectrum  curve  at  energies  less  than  300!  . 
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Another  problem  arises  with  how  large  a  number  the  com¬ 
puter  can  handle.  In  particular,  a  program  execution  error 
occurs  if  the  ratio  £79  in  the  exponent ial  term  in  the 
denominator  of  0(£;9)  is  too  large  (approximately  709  encoun¬ 
tered  in  this  study).  This  condition  occurs  at  all  energies 
when  9«k2T  .  Again,  for  the  forms  of  C( 9)  presented  in  this 
study,  and  at  any  fixed  E  in  the  range  0<£<309t  ,  the 

integrand  of  equation  (20)  adds  a  negligible  contribution  to 
r(£)  at  temperatures  below  about  9  =  £/50 

Taking  into  account  the  above  observations,  the  numer¬ 
ical  approximation  for  equation  (20)  is: 


r»  i 


f( E)  «  I  P(E ;0)C(0)d0  ,  ef  =  £7S0>0 


(23) 


Note  that  the  lower  temperature  limit  changes  with  £  .  The 
approximation  for  equation  (21)  is: 


£(£) 


=  HE1  m  f(f) 


r  rs  i*w 

f(£)d£ 

Js  blah 


*  ^  low  ^  ^  • 


£„,„„  =  309,  <*>  (24) 


For  the  above  two  equations, 
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W'V  -  jp 


ef/9-  l 


(25) 


and  C(6)=di//a«  take  on  the  specific  forms  addressed  in  the 
next  several  chapters. 
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IV. 


Model  X  “  Rad i at i ve  Coo  1 i ng  with  Constant . 


Mater i a  1  Heat  Caoac i t v .  but  No  Rad i at i on  Heat  Capac i t v 

In  this  chapter,  we  look  at  a  description  of  the  model; 
develop  a  specific  form  of  heat  capacity  to  be  used  in  equa¬ 
tion  03);  predict  the  t ime-i ntegrated  spectrum  behavior 
with  a  collection  of  one-temper ature  Planckians;  and  present 
and  discuss  the  results  of  calculating  the  t i me- i ntegrated 
spectrum  for  three  different  initial  temperatures. 

Model  1  Description 

The  first  model  deals  only  with  the  material  particles 
losing  kinetic  energy.  That  kinetic  energy  is  converted  to 
radiation  energy  which  radiates  into  space  from  the  surface 
of  the  plasma  volume.  It  does  so  at  a  rate  that  corresponds 
to  the  release  of  an  equal  amount  of  energy  for  every  degree 
drop  in  temperature.  In  other  words,  the  material  heat 
capacity  is  constant;  it  is  independent  of  temperature.  The 
internal  radiation  density,  on  the  other  hand,  must  be  neg¬ 
ligible  compared  to  the  material  energy  density  to  assume 
there  is  little  or  no  radiation  heat  capacity,  that  is  to 
say,  very  little  loss  of  internal  radiation  energy  per 
degree  change  in  temperature.  In  the  absence  of  internal 
radiation  energy  ((/rC9)"0)  and  internal  energy  to  ionize  to 
any  degree  most  of  the  atoms  <  0  ,  where  "be" 

stands  for  binding  energy),  this  model,  in  general,  repre¬ 
sents  cooling  behavior  at  temperatures  much  lower  than  about 
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0.001  KeV  (on  the  order  of  100,000  K),  the  temperature  at 
which  £/r(9)  and  9)  become  significant  compared  to  the  mate¬ 

rial  energy,  £/„,( 9)  .  This  will  become  evident  in  chapters  VI 
and  VII  where  £/,( 9)  and  are  discussed  relative  to  i/«(9)  , 

respectively.  However,  the  model  is  extended  to  higher  tem¬ 
peratures  to  observe  the  cooling  behavior  of  this  limiting 
case.  This  simple  model  is  used  to  demonstrate  the 
methodology  and,  as  a  limiting  case,  to  show  how  very  non- 
Planckian  the  fluence  spectrum  can  be. 

Heat  Capacity 

The  expression  for  C(8)  ,  or  actually  just  C  since  the 
heat  capacity  in  this  case  is  independent  of  temperature,  is 
derived  from  the  expression  for  the  material  energy  only, 
(/.(«)  in  equation  OS).  Assuming  the  material  particles 
exhibit  a  Maxwel l-Boltzmann  distribution  in  energy,  the 
average  kinetic  energy  of  any  particle  is  (3/ 2)kT  =  (3/ 2)8 

(3:12).  Multiplying  this  by  the  total  number  of  particles, 
and  dividing  by  the  volume  gives  the  material  energy  den¬ 
sity: 

tt.(0)  -  Um(Q)/V  =  Nv( 3/2)9  (26) 

where  uB  is  the  material  energy  density;  and  N v  is  the 
atomic  number  density.  Taking  the  derivative  with  respect 
to  temperature  gives  the  heat  capacity  per  unit  volume: 

Cv  -  CIV  =  dura/d0  =  A/v( 3/2)  (27) 
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which  is  constant,  provided  the  number  density  remains  con¬ 
stant  throughout  the  cooling  process.  (Ionization  is  disre¬ 
garded  until  Model  4).  Replacing  Cv  for  C( 8)  in  equation 
(S3),  it  is  possible  to  bring  Cv  outside  the  integral  sign 
as  a  constant  so  that  equation  < S4 )  becomes: 


*(£) 


cv  f%(£;9)de 


Cv 


f>(E ;  G )  d  9  d£ 


(28) 


Because  Cv  cancels  in  the  above  equation,  it  can  assume  any 
constant  value  without  loss  of  generality.  Therefore,  C,-I 
to  simplify  calculation  of  equation  (38). 


Planckian  Spectra 

Note  that  with  Cr ■ I  ,  the  numerator  of  equation  (38) 
becomes  a  sum  of  normalized  Planckians  at  different  tempera¬ 
tures  over  the  full  range  of  E  .  Figure  3  displays  several 
normalized  Planckian  spectra  at  various  temperatures 
calculated  from  equation  (35).  This  figure  serves  to  show 
the  expected  trend  for  the  curve  of  the  fluence  spectrum. 
Curve  e,  which  is  not  normalized,  is  the  sum  of  curves  a 
through  d,  and  represents  the  approximate  shape  of  the  flu¬ 
ence  spectrum.  It  shows  that  as  the  temperature  continues 
to  drop,  a  greater  percentage  of  the  total  energy  is  carried 
away  by  lower  energy  photons  or,  equivalently,  lower  fre¬ 
quency  photons. 
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Fig.  2.  Planckian  Spectra  a  -  d 
Normalized  to  Unit  Area 


Mathematically  speaking,  the  curve  should  continue  to 
infinity  by  a=0  for  the  fallowing  reason.  If  a=G  in 
equation  <25),  then  the  portion  of  the  equation  equal  to 
i5/Cft4exp(S70)-  i)  remains  constant,  even  as  both  E  and  G 
approach  zero.  However,  G4  in  the  denominator  becomes  the 
dominant  factor  causing  £>  to  approach  infinity  as  it 
approaches  zero.  For  any  fixed  G  greater  than  zero,  and  as 
E  approaches  zero,  the  fi"3  term  in  the  numerator  decreases 
faster  than  the  denominator,  and  is  the  reason  why  all  of 
the  above  spectra  diminish  to  zero  by  E=0  .  Adding  more 

Planckians  to  Figure  2  would  enlarge  curve  e,  but  the  curve 
should  retain  its  approximate  shape.  Lower  temperature 
Planckians  added  successively  to  the  left  of  curve  d  would 
add  diminishing  areas  to  the  total  area  underneath  curve  e. 


19 


Thus,  when  curve  e  is  normalized  to  unit  area,  the  peak 
should  reach  some  finite  limit  at  the  left  edge  of  the 
graph . 

Time-Integrated  Spectra 

Time-integrated  spectra  with  initial  temperatures  of 
0.02,  1,  and  10  KeV  are  calculated  from  equation  (28)  using 

the  computer  program  listed  in  Appendix  A. 

Figure  3  shows  the  time-integrated  spectrum  with  an  ini¬ 
tial  temperature  of  1  KeV . 


Fig.  3.  Model  1  Time-Integrated  Spectrum 

Normalized  to  Unit  Area;  Q ,  =  l  KeV 

Indeed,  the  spectrum  curve  has  roughly  the  same  shape  as 
curve  e  in  Figure  2  from  the  peak  to  the  upper  end  of  the 
spectrum,  and  the  peak  approaches  a  limiting  value  of  about 
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0.37  relative  units  of  energy  density.  Upon  closer  examina¬ 


tion  at  the  left  edge  of  the  graph  shown  in  Figure  4,  the 
spectrum  curve  rolls  over  at  about  £*0.01  to  drop  to  just 
less  than  0.05  near  £  =  0 


0  0.02  0.04  0.06  0.08  0.1  0.12  0.14  0.16  0.18  0.2 

Photon  Energy  (KeV) 


Fig.  4.  Model  1  Time-Integrated  Spectrum 
(Magnified  View)  Normalized  to 
Unit  Area;  9t  =  1  KeV 

Actually,  £  =  0.001  was  the  value  for  the  lower  integration 
limit  in  equation  (28)  for  this  example.  The  curve  would  in 
fact  continue  to  approach  zero  at  the  left  end  as  E 
approaches  zero.  The  use  of  £=0/1000  instead  of  another 
value  even  closer  to  zero  is  justified  because  the  approxi¬ 
mately  triangular-shaped  area  to  the  left  of  the  curve  from 
the  ordinate  value  of  0.37  down  to  zero  would  be  a 
negligible  addition  to  the  area  under  the  curve  <1/2  X  0.37 
X  0.01)  =  0.00185  or  about  0.2'/.  compared  to  1).  The  small 
peak  at  the  left  end  in  Figure  4  is  an  artifact  of  the 
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calculation  that  becomes  more  narrow  and  displaced  to  the 
left  as  the  temperature  mesh,  A0,  is  made  smaller,  so  in 
essence,  it  may  be  disregarded. 

Spectra  for  initial  temperatures  of  0.02  and  10  KeV  are 
shown,  respectively,  in  Figures  5  and  6.  Comparison  of  all 
three  time-integrated  spectra  reveals,  in  general,  that  the 
spectrum  shifts  toward  lower  energies  for  lower  initial  tern 
peratures  and  toward  higher  energies  for  higher  initial  tern 
peratures,  the  same  as  for  single-temperature  Planckians. 
For  instance,  considering  the  areas  under  the  curves, 
practically  100*/.  of  the  total  energy  is  carried  away  from 
the  plasma  by  photons  in  the  range  0  to  1  KeV  when  9,-0. 02 
KeV,  about  33V.  of  the  energy  when  9,-1  KeV,  and  only  3*/. 
when  9,-10  KeV.  Another  feature  of  the  time-integrated 
spectra  similar  to  Planckian  spectra  is  the  fact  that  near  1 
all  of  the  internal  energy  is  radiated  away  by  photons  hav¬ 
ing  energies  less  than  100,  KeV. 


[•It 


Fig.  5.  Model  1  Timp- Integrated  Spectrum 

Normalized  to  Unit  Area;  8,- 10  KeV 
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Photon  Energy  (KeV) 

Fig.  6.  Model  1  Time-Integrated  Spectrum 

Normalized  to  Unit  Area;  9|  =  0.02  KeV 

Significant  differences  between  the  two  types  of  spectra 
are  the  fraction  of  energy  radiated  at  the  lower  end  of  the 
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spectrum  and  the  most  probable  energy  for  photons.  For  any 
Planckian,  approximately  1/30  of  the  total  plasma  energy  is 
radiated  at  energies  below  its  temperature,  9  KeV,  whereas 
for  any  Model  1  time-integrated  spectrum,  about  ten  times 
the  Planckian  fraction  or  one  third  of  the  total  energy  is 
radiated  at  energies  below  its  initial  temperature,  9,  KeV. 
Over  time,  the  most  probable  energy  for  photons  of  a  con¬ 
stant  temperature  Planckian  emitter  is  £  =  2.829  ,  corre¬ 

sponding  to  the  energy  location  of  the  peak.  For  the 
cooling  Planckian  emitter,  the  most  probable  energy  for 
photons  shifts  toward  lower  energies  over  time,  the  time- 
integrated  spectrum  showing  it  near  E= 0 

Summary 

Model  1  represents  an  upper  bounding  case  for  energy 
density  at  the  low  end  of  the  spectrum  for  the  cooling 
Planckian  emitter.  With  all  of  its  energy  in  the  form  of 
material  kinetic  energy  and  the  heat  capacity  remaining  con¬ 
stant  throughout  the  cooling  process,  there  is  an  upper 
limit  to  the  fraction  of  total  energy  radiated  at  the  lower 
end  of  the  spectrum.  It  turns  out  that  density  does  not 
affect  the  normalized  spectrum,  as  long  as  there  exits  some 
amount  of  material.  If  material  energy  constitutes  a  sig¬ 
nificant  portion  of  the  total  energy  of  a  cooling  body,  its 
most  influential  effect  on  the  t i me- i ntegr a ted  spectrum 
would  be  seen  at  lower  temperatures  at  the  lower  end  of  the 
spectrum . 
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In  the  next  chapter  we  look  at  a  model  of  the  cooling 
plasma  containing  only  radiation  energy  to  observe  its 
effects  upon  the  fluence  spectrum. 
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V.  Model  2  -  Rad i at i ve  Coo  1 i no  with  Temperature  Cons i stent . 

Rad i at i on  Hea t  Capac i tv .  but  No  Mater i a  1  Heat  Capac i t v 

We  looked  at  the  effect  of  including  only  material 
energy  in  the  cooling  emitter  in  the  last  chapter.  Now  we 
take  a  look  at  the  effect  of  including  only  radiation 
energy . 

Model  2  Description 

The  second  model  deals  with  the  plasma  cooling  through 
the  loss  of  radiation  energy  only.  Either  the  material  par¬ 
ticles  are  at  rest,  with  their  kinetic  energy  being  close  to 
zero  (but  this  would  assume  the  body  has  already  cooled  to  a 
temperature  near  0  KeV),  or  else  the  material  energy  density 
is  negligible  compared  to  the  radiation  energy  density  at 
higher  temperatures ,  which  is  more  likely  the  case.  At  very 
high  temperatures,  well  above  the  point  where  the  atoms  com¬ 
pletely  ionize,  a  drop  in  temperature  to  a  point  still  above 
complete  ionization  would  result  in  no  change  in  heat 
capacity  due  solely  to  ionization.  So  this  model  seems 
appropriate  at  very  high  temper atures .  Nevertheless,  the 
model  is  extended  to  low  temperatures  to  observe  the  cooling 
behavior  there  also,  and  binding  energy  is  not  taken  into 
account . 
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Heat  Capacity 


The  heat  capacity  in  this  approximation  is  dependent  on 
temperature  because  of  the  high  degree  of  dependence  that 
£/r(9)  has  on  temperature.  The  form  for  radiation  energy  den¬ 
sity  is: 

ur(9)  -  Ur(Q)/V  =  a94  (29) 


where  a  is  the  radiation  constant  equal  to  3 . 56E+3B  in  units 
of  l/(KeV-m)3  (3:32).  The  heat  capacity  is  then: 

Cw(0)  -  C(0)/K  =  du,(0)/d0  =  4a03  (30) 


Replacing  C( 9)  in  equation  (33)  with  C»( 9)  and  bringing  the 
constant,  4a  ,  outside  the  integral  in  the  numerator  and 
denominator  of  equation  (34),  it  is  seen  that  the  constants 
cancel  to  give: 


/-«< 

/  A(f;e)0: 


d9 


/•futon  rsi 

£(£  ;0)03d9  d  E 

J  S  l#w 


(31) 


Note  that  Cv(9)  could  just  as  well  be  9s  without  loss  of 
generality,  which  is  the  form  of  heat  capacity  used  to  cal¬ 
culate  the  fluence  spectrum  for  Model  2. 

Because  the  heat  capacity  has  a  0*  dependence,  more 
energy  is  released  from  the  plasma  volume  per  degree  drop  in 
temperature  at  higher  temperatures  than  at  lower  ones.  Pho¬ 
tons  of  higher  energy  should  carry  away  more  of  the  total 
energy  than  in  the  case  of  constant  material  heat  capacity. 
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Correspondingly,  the  peak  of  the  spectrum  should,  in  gen¬ 
eral,  shift  toward  higher  energies,  but  the  overall  spectrum 
behavior  must  be  determined  from  the  heat  capacity  in 
conjunction  with  the  function  £;0)  . 

Time-Integrated  Spectra 

Predicting  the  shape  of  the  spectrum  can  be  done  by 
examining  the  integrand  in  the  numerator  of  equation  (31) 
with  £(£;9)  defined  in  equation  ( 25 )  : 

$(*:&)  e3  = 

ft4 

In  the  limiting  case  at  the  low  end  of  the  spectrum,  with 
£“9  and  £  approaching  zero,  the  entire  expression  is  con¬ 
stant  except  for  a  remaining  E2  term  in  the  numerator. 
Therefore,  as  E2  pproaches  zero,  the  integrand  approaches 
zero.  At  the  'igh  end  of  the  spectrum  with  9<9,<®  ,  and  E 

approaching  infinity,  the  exponential  in  the  denominator 
approaches  infinity  faster  than  the  numerator;  hence,  the 
integrand  again  approaches  zero.  In  between  the  two  limits, 
a  peak  should  occur  at  the  same  most  probable  energy  as  a 
Planckian  with  fixed  temperature  9  (£  =  2.828),  and  the  entire 
curve  is  scaled  in  height  by  the  parameter  83  ,  as  seen  in 
the  left  hand  side  of  equation  (32). 

Figure  7  shows  the  effect  of  adding  several  of  the  func¬ 
tions  given  by  equation  (32)  for  a  few  different  tempera¬ 
tures.  Curve  d  simulates  approx imately  what  the  fluence 
spectrum  might  look  like  for  a  cooling  body  with  a  starting 


—  -* 
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temperature  of  1  KeV .  Notice  that  the  peak  has  shifted  to  a 
higher  energy  as  discussed  above,  but  is  left  of  the  peak 
for  a  Planckian  with  temperature  8,  .  The  second  condition 
follows  from  the  fact  that  as  the  body  cools  the  energy 
distribution  shifts  to  lower  energies.  Photons  are  most 
likely  to  have  energies  that  are  nearly  equivalent  to  the 
body  temperature. 


Fig.  7.  Planckian  Spectra  a  -  c  Normalized 
to  Unit  Area  and  Multiplied  by  9’ 


Furthermore,  the  second  condition  means  there  is  an  upper 
limit  to  the  fraction  of  total  energy  that  can  be  radiated 
at  the  high  end  of  the  spectrum.  For  example,  a  Planckian 
with  temperature  0,  KeV  radiates  a  maximum  of  about  95.67.  (or 
29/30)  of  the  total  energy  at  photon  energies  higher  than  0, 
KeV.  This  upper  limit  should  hold  true  for  any  time- 
integrated  spectrum. 
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Fig.  9.  Model  2  Time-Integrated  Spectrum 

Normalized  to  Unit  Area;  8,-0. 02  KeV 


Fig.  10.  Comparison  of  Time-Integrated  Spectrum 
F=1  KeV  with  Plankian  Spectra 


All  three  t ime- integrated  spectra  have  the  same  shape  on 
a  relative  scale,  but  are  located  at  different  points  if 
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placed  on  the  same  scale.  The  same  trend  occurs  as  in  Fig¬ 
ure  1  where  the  hotter  spectra  have  more  of  their  energy 
density  at  higher  photon  energies.  Their  peaks  ,ire  ■smaller 
to  accommodate  the  same  area  under  the  curves,  and  shifted 
to  the  right,  to  reflect  the  most  probable  energies  being 
within  a  factor  of  three  from  the  starting  temperatures. 

Just  the  opposite  is  true  of  the  colder  spectra;  the  peaks 
grow  taller  and  shift  to  the  left  for  the  same  reasons. 

Comparison  with  Plankian  Spectra 

Consider  the  t ime- i ntegr ated  spectrum  in  Figure  10  with 
a  starting  temperature  of  1  KeV.  The  shape  of  the  spectrum 
is  similar  to  a  Planckian,  but  as  anticipated,  the  peak  is 
displaced  to  the  left  of  that  for  a  1  KeV  Planckian  (the 
Planckian  shown  is  one  having  9=0,7345  KeV).  The  peak  for 
the  time-integrated  spectrum  is  located  at  E  = Z.079|  com¬ 

pared  to  E- Z,5Z9|  for  the  corresponding  Planckian  peak,  and 
the  cooling  body  radiates  about  92%  of  its  total  energy  at 
energies  above  E- 9(  compared  to  96'/.  for  the  Planckian 
emitter  always  at  a  constant  temperature. 

Figure  10  shows  the  difference  in  shape  between  both 
types  of  spectra  having  both  their  peaks  at  the  same  energy, 
5  .  It  can  be  seen  that  one  Planckian  roughly  fits  the 
t i me- i ntegr a ted  spectrum,  whereas  a  combination  of  two 
Planckians  with  different  temperatures  practically  matches 
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the  spectrum.  Appendix  B  contains  a  computer  program  to  do 
the  curve  fitting  using  two  Planckian  basis  functions.  The 
equations  for  the  fit  are  given  by: 

Pi  CF)  =  FVCexpCF/O,)-!]  (33) 

Fa(F)  =  F3/[exp(£/0a)- I]  (34) 

*m(*>  =  aiPl(iE)  +  a2P2(E)  (35) 

The  measure  of  the  fit  is  given  by  the  sum  of  the  square  of 
the  errors: 

SSE  =  £>(£)  -  $,lt(£))2  (36) 

t 

where  a(£)  is  the  time-integrated  spectrum  value  from  equa¬ 
tion  (24),  and  SntCfi’)  is  the  curve-fit  value  given  by  equa¬ 
tion  (35).  From  this  result,  the  cooling  body  can  be 
thought  of  as  radiating  at  two  temperatures ,  both  below  the 
initial  temperature,  instead  of  many  temperatures  over  the 
course  of  time. 

Summary 

Model  2  represents  an  extreme  bound  in  the  opposite 
direction  to  Model  1.  It  shows  there  is  an  upper  limit  to 
the  fraction  of  total  energy  that  can  be  radiated  at  the 
high  end  of  the  spectrum.  If  the  cooling  body  contains  a 
significant  portion  of  radiation  energy,  its  effect  would  be 
seen  at  higher  temperatures  and  energies  than  the  case  for 
constant  material  heat  capacity,  and  in  fact  should  be  the 
dominant  factor  at  the  high  temperatures  and  energies  under 
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VI .  Mode  1  3  -  Rad i at i ve  Coo  lino  wi th  Constant .  Mater i a  I 
Heat  Capac i tv  and  Temperature  Cons i stent .  Rad i at  ion  Heat 

Caoac i tv 

This  chapter  combines  the  effects  of  both  material  and 
radiation  energy  in  the  cooling  body.  In  the  previous  two 
models,  temperature  is  the  only  parameter  that  influences 
the  energy  distribution,  but  now  density  becomes  important, 
as  we  shall  see. 

Model  3  Description 

The  third  model  for  the  cooling  body  takes  into  account 
the  loss  of  kinetic  energy  from  the  material  particles  and 
the  loss  of  radiation  energy  or  those  photons  that  escape 
any  further  interaction  withir,  the  plasma.  This  model  could 
be  valid  at  low  temperatures  and  at  very  high  temperatures, 
but  excludes  the  range  of  intermediate  temperatures  where 
the  ionization  state  of  the  atoms  is  changing  with  tempera¬ 
ture.  Binding  energy  is  not  included,  but  the  intermediate 
range  of  temperatures  is,  so  that,  as  before,  cooling 
behavior  can  be  observed  here  too.  The  densities  considered 
may  be  somewhat  high  for  plasmas  in  fusion  experiments;  how¬ 
ever,  they  are  chosen  to  illustrate  noticeable  effects  for 
spectra  with  initial  temperatures  of  0.03  and  1  KeV. 


35 


Heat  Capacity 


Both  the  formula  for  material  energy  density  from  Model 
1  and  the  formula  for  radiation  energy  density  from  Model  2 
comprise  the  total  internal  energy  density  for  Model  3: 

u(0)  -  i/(0)/ V  =  um(0)  +  u,(0)  (37) 

=  Nr(3/2)9  +  ct04  (38) 

Notice  here  why  the  effect  due  to  material  energy  is  domi¬ 
nant  at  low  temperatures,  and  the  effect  due  to  radiation 
energy  is  dominant  at  high  temperatures.  Both  A/„  and  a  are 
about  the  same  order  of  magnitude  for  nominal  material  den¬ 
sities,  but  04  is  much  less  than  0  for  8<«1  KeV,  and  it  is 
much  greater  than  0  for  0^1  KeV.  The  two  terms  are  equal 
at  approximately  0-0.56  KeV  (6.5  million  K>  when  JVr-lEf28 
atoms/m53 . 

The  derivative  of  equation  (38)  with  respect  to  tempera¬ 
ture  is  the  expression  for  heat  capacity  for  Model  3: 

C¥(0)  -  C(0)/K  -  d  U(0)/d9  (39) 

=  Nv(  3/2)  +  4a03  (40) 


Time-Integrated  Spectra 

Substituting  C'f( 8)  for  C(8)  in  equation  (24)  results  in 
the  equation  to  compute  the  time-integrated  spectrum  for 
this  model : 


&<£)  - 


I  P(£;B)C'(B)  d0 

Jit 


rjston  rst 

p(£ ;0)Cv(0)d0  d E 

J  t  t»w  Jit 


(41) 
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All  the  spectra  for  Model  3  are  calculated  from  equation 
(41)  using  the  computer  program  in  Appendix  A,  with  the  sub¬ 
routine  to  compute  heat  capacity  changed  to  reflect  CT(9)  in 
equat ion  ( 40 )  . 

The  predicted  cooling  behavior  depends  on  the  tempera¬ 
ture.  At  high  initial  temperatures,  most  of  the  energy  loss 
is  due  to  the  change  in  internal  radiation  energy. 

Therefore,  the  heat  capacity  has  an  almost  pure  93  depen¬ 
dence  like  that  of  Model  2,  so  the  spectrum  should  resemble 
a  Model  2  spectrum  with  the  same  high  initial  temperature. 
The  peak,  however,  should  be  displaced  more  to  the  left 
because  as  the  plasma  cools,  more  energy  is  lost  from  mate¬ 
rial  than  would  be  the  case  for  radiation  only.  Such  is  the 
case  for  the  10  KeV  spectrum  of  a  plasma  with  density  equal 
to  1E+27  atoms/m3;  however,  the  shift  is  too  small  to  show 
any  appreciable^jlifference  with  the  Model  2  spectrum  in  Fig¬ 
ure  8.  For  low  initial  temperatures ,  most  of  the  energy 
loss  is  due  to  the  change  in  material  energy.  The  heat 
capacity  is  nearly  conut. mt  and  becomes  more  so  as  the 
plasma  cools.  In  this  case,  the  spectrum  should  resemhln  ,j 
Model  1  spectrum  having  the  same  initial  temperature,  but 
with  the  peak  shifted  right  because  of  the  energy  losses 
from  internal  radiation  at  temperatures  closer  to  the  ini¬ 
tial  temperature.  The  0.02  KeV  spectrum  behaves  thus,  but 
again  has  a  shift  too  small  to  show;  it  is  almost  identical 
to  the  spectrum  in  Figure  9  in  Chapter  V. 
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Cooling  behavior  from  intermediate  starting  temperatures 
should  show  noticeable  effects  due  to  both  material  and 
radiation  energy.  The  peak  should  fall  between  those  of  the 
pure  material  case  and  the  pure  radiation  case  for  the  same 
starting  temperature.  Figure  11  gives  an  example  of  this 
behavior,  with  the  Model  3  spectrum  calculated  for  a  plasma 
density  of  1E+29  atoms/m3. 


Fig.  11.  Model  1,  2,  and  3  Ti me- I ntegr ated  Spectra 
Normalized  to  Unit  Area;  Qj  =  i  KeV 

Note  that  the  left  end  of  the  Model  3  spectrum  approaches  a 
value  greater  than  zero  because  of  the  material  contribution 
to  the  loss  of  energy.  This  becomes  more  pronounced  for 
colder  9,  ,  while  for  hotter  9,  the  value  approaches  zero. 

The  predicted  cooling  behavior  depends  on  density  as 
well.  Higher  densities  mean  more  material  energy  which 
forces  the  heat  capacity  to  become  more  constant  in  nature 
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beginning  at  higher  temperatures.  The  effect  of  increasing 
the  density  causes  the  lower  end  of  the  spectrum  to  approach 


its  peak  limiting  value  for  pure  material  energy,  as  shown 
in  Figure  12.  A  density  of  1E+29  (not  shown)  versus  1E+27 
atoms/m3  resulted  in  no  noticeable  change  in  the  spectrum. 
See  Figure  6  in  Chapter  IV  for  the  limiting  case. 


Fig.  12.  Density  Effect  on  Model  3  T i me- I ntegr a ted 
Spectra  Normalized  to  Unit  Area; 
e,  =  Q.02  KeV 

Lower  densities  mean  less  material  energy,  and  the 
nature  of  the  heat  capacity  becomes  more  like  the  case  for 
pure  radiation.  Figure  13  is  a  better  example  of  what  low¬ 
ering  the  density  does;  the  spectrum  becomes  more  like  the 
pure  radiation  case  shown  in  Figure  11  ,  or  Figure  10  in 


Chapter  V. 


Fig.  13.  Density  Effect  on  Model  3  Time-Integrated 
Spectra  Normalized  to  Unit  Area; 

9t=  1  KeV 


Comparison  with  Planckian  Spectra 

Chapter  V  showed  that  two  Planckian  basis  functions  fit 
the  fluence  spectra  of  Model  5  much  better  than  a  single 
basis  function.  This  is  generally  true  for  any  spectrum 
having  a  shape  closely  resembling  that  for  the  pure  radi¬ 
ation  case.  However,  the  curve-fit  becomes  worse  for  low 
temperature  or  high  density  Model  3  spectra  as  they  become 
similar  in  shape  to  Model  1  spectra.  Figures  14  and  15 
support  these  conclusions  for  different  densities. 
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Photon  Energy  (KeV) 


.  14.  Comparison  of  Two-Pl anck i an  Curve-Fit 

with  Model  3  Spectrum;  My~lE+27 
particles/m3 


0123456789  10 

Photon  Energy  (KeV) 


Fig.  15.  Comparison  of  Two-Planck ian  Curve-Fit 
with  Model  3  Spectrum;  Nr - l E+29 
part ic les/m3 
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Summary 


Model  3  represents  a  more  realistic  theory  of  cooling 
behavior  between  the  extreme  bounds  of  Model  1  and  2.  The 
fluence  spectrum  generated  by  Model  3  combines  effects  of 
both  material  and  radiation  energy  losses  which  place  the 
most  probable  energy  somewhere  between  the  two  extremes, 
depending  on  density  and  initial  temperature  of  the  plasma. 
In  the  case  of  high  densities  or  low  temperatures,  Model  3 
degenerates  into  the  same  behavior  as  Model  1;  the  heat 
capacity  is  nearly  constant.  At  low  densities  or  high  tem¬ 
peratures,  it  becomes  like  Model  2;  the  heat  capacity 
behaves  as  9a  .  Although  more  realistic,  this  model  is 
simplistic  and  incomplete  without  taking  into  account  the 
release  of  energy  as  the  ionization  state  falls  for  the  ions 
in  the  cooling  plasma.  This  aspect  is  covered  in  the  next 
chapter . 
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VI  I 


Node  1  4  -  Rad i at i ve  Coo  1 i no  with  Temperature  Consi s- 


tent  Heat  Capac i tv 

Model  4  concentrates  on  how  material  composition  affects 
the  fluence  spectrum  of  the  cooling  Planckian  emitter.  The 
form  of  material  energy  presented  in  Model  1  is  changed  to 
reflect  the  variable  number  of  free  electrons  present  in  the 
plasma  as  the  temperature  falls  and  they  bind  to  the  ions. 
Binding  energy  is  accounted  for  in  yet  another  internal 
energy  term  to  calculate  the  heat  capacity.  The  results 
include  variation  in  initial  temperature  and  type  of  single 
element  material,  but  the  general  effect  due  to  variation  in 
initial  density  is  left  to  the  previous  chapter. 

Model  4  Description 

The  fourth  model  describes  the  cooling  plasma  through 
the  loss  of  three  types  of  internal  energy:  material 
energy,  radiation  energy,  and  binding  energy,  all  functions 
of  temperature.  As  the  temperature  falls,  free  electrons  in 
the  plasma  lose  kinetic  energy  and  are  captured  through 
electrostatic  attraction  by  the  ions,  which  eventually  fill 
all  electron  shells  to  become  neutral  atoms.  The  shell 
structure  is  governed  by: 

nt,  =  2  jz  .  7=1  .2....  (42) 

where  is  the  total  number  of  electrons  allowed  in  princi¬ 
pal  quantum  state  j  .  Electrons  are  assumed  to  fill  inner 
vacancies  first  and  move  progressively  away  from  the  nucleus 
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to  fill  outer  vacancies  last.  When  an  electron  is  captured, 
it  proceeds  directly  to  the  next  unfilled  vacancy,  emitting 
a  photon  of  energy  equal  to  the  absolute  difference  between 
the  energy  level  of  the  vacancy  just  filled  and  the  contin¬ 
uum.  No  allowance  is  made  for  intermediate  excited  states. 
The  entire  population  of  ions  is  considered  to  be  at  the 
same  ionization  state;  however,  fractional  numbers  of  bound 
electrons  are  permitted  to  account  for  a  distribution  of 
ionization  states  encountered  in  real  physical  situations. 
This  model  applies  primarily  for  all  temperatures  above  the 
point  where  the  material  vaporizes.  But  as  long  as  the 
initial  temperature  is  hot  enough  to  include  some  degree  of 
ionization,  then  energy  losses  attributable  to  latent  heats 
of  vaporization  and  fusion  are  negligible  compared  to  bind¬ 
ing  energy,  and  the  model  is  valid  at  all  temperatures. 

There  is  a  restriction,  however,  on  the  upper  limit  for 
density,  depending  on  type  of  material,  because  Model  4  uses 
a  form  of  the  Saha  equation  that  assumes  the  free  electrons 
are  nondegenerate. 
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r  j  _  ,  (Z;Y  .  n-*Y  18  Cr«), 

Ej  -  HjJ  2*  It  -  in- 


1.2....  (43) 


where  /«  is  the  ionization  potential  for  an  electron  in  the 
ground  state  of  a  hydrogen  atom;  Z*  is  the  effective  nuclear 
charge  as  seen  by  th«-  electron  in  the  j til  energy  state;  n, 
is  the  number  of  free  electrons  per  atom;  s  is  the  elec¬ 
tronic  charge  of  an  electron;  R  is  the  radius  of  the  spheri¬ 
cal  region  surrounding  the  ion  which  contains  zero  net 
charge,  i.e.,  n,  free  electrons;  and  (rz),  is  the  mean  square 
radius  of  the  electron  orbit.  The  left  hand  term  gives  the 
energy  levels  for  a  hydrogen-like  atom  with  atomic  number 
Z*,  while  the  right  hand  term  corrects  for  interactions 
between  bound  and  free  electrons,  free  electrons  with  the 
ion,  and  free  electrons  with  free  electrons,  and  adjusts  the 
zero  energy  level  so  that  all  free  electrons  within  the  ion 
sphere  maintain  positive  energy.  (Bound  electrons  have  neg¬ 
ative  energy).  For  convenience,  equation  (43)  computes  E , 
in  units  of  eV  instead  of  Rydbergs  if  !„  =  13.61  eV  and  the 
right  hand  term  is  multiplied  by  l/[4Tvf«Cl  .6C2E~  19  J/SV)]  , 
with  g0  ,  the  vacuum  permi ttivi ty ,  in  F/m  and  R  in  meters. 
The  equation  for  effective  nuclear  charge  is  (3:168): 


r  i  \ 

z  -  Z'1***,  -  n,{1  - 


(44) 


where  Z  is  the  element  atomic  number;  nt  and  n-,  are  the 
number  of  electrons  in  the  corresponding  shells  scrpRiinq 
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the  nuclear  charge  from  the  electron  in  question;  and  at)  and 
o„  are  Slater  screening  constants  (3:169).  To  find  the  mean 
square  radius  (3:170): 


?Ul[l  + 

(z;)2U 


(45) 


where  a0 
« 

meters . 
given  by 


is  the  first  Bohr  radius  for 
Finally,  the  ion's  sphere  of 
the  relationship  (3:170): 


the  hydrogen  atom  in 
influence  has  a  radius 


(46) 


where  Nv  is  the  atomic  number  density. 

From  equation  (43)  and  (46)  it  is  seen  that  density  is 
one  factor  that  determines  the  binding  energy,  , 

required  to  remove  the  electron  from  the  ion  to  become  a 
free  electron.  A  plot  of  binding  energy,  ~Ej  ,  versus  den¬ 
sity,  ,V*  ,  in  Figure  16  for  the  outer  shell  of  electrons  in 
carbon  reveals  that  only  densities  less  than  1E+E8  atoms/m3 
satisfy  the  requirement  that  it  takes  more  energy  to  remove 
the  next  electron  than  the  previous  one.  Above  this  limit, 
the  electron  density  is  so  high  that  the  electrons  become 
degenerate;  that  is,  more  than  one  electron  occupies  the 
same  energy  level.  This  upper  density  limit  differs  with 
atomic  number  Z  ,  but  for  the  four  elements  chosen,  carbon, 
aluminum,  copper,  and  plutonium,  Nv  =  1E+27  atoms/m3  is 
satisfactory  for  comparison  of  their  spectra. 
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Fig.  16.  Binding  Energy  versus  Atomic  Number 
Density  for  Carbon 


The  binding  energy  required  to  remove  the  outermost 
electron  from  an  ion  or  neutral  atom  can  be  related  to  tem¬ 


perature  by  the  following  expression,  taken  from  reference 
(4),  which  is  derived  from  the  Saha  equation  (3:164)  for  a 
plasma  containing  a  single  Z  material: 


2  f  2nmtQ 


'111  - 

- 


N  t  =  ntN , 


/  =  1.2....  (47) 


(48) 


n,  =  Z  -  nt  -  ^2(y-l)i 


(49) 


where  m,=  5.685F- 15  KeV-ss/ma  is  the  mass  of  an  electron;  N, 
is  the  electron  density  in  1/m3  ;  ft=4.l36fT-18  KeV-s  is 

Planck's  constant;  and  Et  and  0  are  in  KeV.  If  N *  and  n( 
are  given,  then  E,  and  9  can  be  calculated  from  equations 
(43)  -  (49).  0  is  the  temperature  that  corresponds  to  nr 
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free  electrons  per  atom.  -Et  is  the  energy  needed  to  free 
the  next  electron  to  be  added  to  the  ft,  free  electrons.  It 
i  not  included  in  the  total  binding  energy  required  to 
remove  ft,  electrons,  that  is,  the  sum  of  the  binding  ener¬ 
gies  to  remove  each  of  the  individual  ft,  electrons. 

So  at  temperature  9  ,  there  are  N»ft,  free  electrons  per 
unit  volume  in  the  plasma,  and  the  total  binding  energy  is: 

U..C9)  -  (/..//  -  -W,£s(Cn ,(8»  (SO) 

where  the  sum  starts  with  0  KeV  for  ft,-G  and  adds  the 
binding  energy  for  each  additional  electron  removed  from  the 
bound  system.  The  function  ft, (9)  is  linearly  interpolated 
from  a  table  of  values  for  ft,  versus  0  calculated  from  equa¬ 
tions  (43)  -  (49)  for  each  element.  A  graph  of  ft,  versus  9 
in  Figure  17  for  copper,  that  is  typical  of  other  elements, 
shows  that  equation  (47)  breaks  down  for  transitions  between 
electron  shells  when  tij  approaches  zero. 

Because  it  is  physically  impossible  to  have  more  than 
one  value  of  ft,  electrons  free  from  the  atom  at  one  tempera¬ 
ture  (which  corresponds  to  one  instant  in  time),  the  graph 
is  smoothed  as  shown  in  Figure  18  by  removing  some  of  the 
points.  In  general,  all  elements  have  this  "S"-shaped 
curve,  and  full  ionization  occurs  at  increasing  temperatures 
for  i  .creasing  2. 
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30 


Fig.  17.  Linear-Log  Plot  of  Free  Electrons  per 
Atom  versus  Temperature  for  Copper 


Fig.  18.  Smoothed  Linear-Log  Plot  of  Free  Electrons 
per  Atom  versus  Temperature  for  Copper 


The  function  of  total  binding  energy  versus  number  of 
free  electrons  per  atom  is  linearly  interpolated  from  a 
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table  also.  Examples  of  both  tables  are  given  in  Appendix  A 
for  copper.  A  typical  graph  of  total  binding  energy  versus 
temperature  is  shown  in  Figure  19.  All  elements  have  this 
irregular,  "S"-shaped  curve,  with  the  abrupt  changes  in 
slope  caused  by  electron  shell  transitions  and  significant 
changes  in  slope  from  the  corresponding  free  electron  versus 
temperature  curves.  In  general,  total  binding  energy 
increases  with  temperature  and  increasing  2. 


Fig.  19.  Linear-Log  Plot  of  Binding  Energy  per 
Atom  versus  Temperature  for  Copper 


Appendix  C  has  two  computer  programs  to  compute:  1)  the 
binding  energy  for  a  particular  electron,  and  S)  the  corre¬ 
sponding  temperature,  from  equations  (43)  -  (49). 


Material  Energy 

The  material  energy  present  in  the  plasma  at  any  given 


time  is  dependent  on  the  temperature  and  the  number  of  free 


particles  which  is  decreasing  because  electrons  bind  to  the 
ions  as  the  temperature  falls.  Therefore,  the  material 
energy  term  is  changed  to  reflect,  in  addition  to  the  ions, 
the  number  of  free  electrons  present  as  a  function  of  tem¬ 
perature,  with  each  particle  having  energy  equal  to  (3/2)9  : 

ura  -  W¥(  1  +  /if(0))(  3/2)9  (51) 


Heat  Capacity 

The  heat  capacity  is  calculated  from  the  expression  for 
internal  energy  density  which  now  has  the  form: 

u(9)  =  N  v(  I  +  ft*  (6))  (3/2)0  +  a04  +  Wv£[-£/(/i,(e))X52) 


Because  nf(5)  is  different  for  every  element  and  somewhat 
difficult  to  express  as  an  analytical  function,  the  deriv¬ 
ative  du/d0  is  approximated  by  the  change  in  internal  energy 
over  a  small  change  in  temperature  A0  as  given  by: 


C  -  —  = 
A  9 


I 


u  9+ -A9  -  u  9- -A9 


I 


/  A  0 


(53) 


To  get  a  feel  for  the  magnitude  of  the  change  in  binding 
energy  compared  to  the  change  in  material  energy  per  temper¬ 
ature  change  of  1  KeV  consider  the  graphs  for  copper  again 
in  Figures  18  and  19.  In  Figure  18,  the  number  of  free 
electrons  at  a  temperature  of  1  KeV  is  about  28,  and  at  0.03 
KeV  it  is  about  7.  The  change  in  material  energy  per  KeV 
per  atom,  then,  is: 
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=  (1  +  28)|(1)  -  (1  +  7)|(0.03)  /  (1-0.03) 

=  44  (KeV/KeV-atom) 

In  Figure  19,  the  binding  energy  at  0=1  KeV  is  about  27  KeV, 
and  at  9-0.03  KeV  it  is  1  KeV.  Hence,  the  change  in  binding 
energy  per  KeV  per  atom  is: 

=  (27-1)  /  (1-0.03)  =  27  (KeV/KeV-atom) 

iV  v 

The  binding  energy  contributes  approximately  60*/.  as  much 
change  in  energy  to  the  heat  capacity  as  does  the  material 
energy.  For  the  other  three  elements,  the  ratio  is  approxi¬ 
mately  in  the  range  of  30*/.  -  110*/.  .  Thus,  the  change  in 
binding  energy  is  the  same  order  of  magnitude  as  the  change 
in  material  energy. 

Time-Integrated  Spectra 

The  t i me-integrated  spectra  for  Model  4  are  calculated 
from  equation  (52)  and  (53)  substituted  into  equation  (41) 
from  Chapter  VI.  Appendix  A  has  the  computer  program  to  do 
this  with  the  appropriate  change  to  Cv  in  the  subroutine  to 
compute  heat  capacity. 

The  contribution  to  heat  capacity  from  material  energy 
in  Model  4  is  an  order  of  magnitude  or  greater,  depending  on 
2,  than  for  Model  3,  at  temperatures  above  0.1  KeV  when  the 
atomic  number  density  ,V,  is  the  same.  Add  to  this  roughly 
the  same  contribution  from  binding  energy,  and  the  influence 
on  the  t i me- i n tegr a ted  spectrum  due  to  the  combined  change 
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in  material  and  binding  energy  is  much  greater  than  Model  3 


The  combined  effect  is  most  significant  in  the  temperature 
range  where  most  of  the  electron  binding  occurs,  generally 
between  0.01  and  1  KeV  for  elements  from  carbon  to  pluto¬ 
nium.  Spectra  having  initial  temperatures  of  1  KeV  should 
show  significant  increases  in  energy  density  near  the  low 
end  of  the  spectrum  with  a  corresponding  shift  in  the  main 
peaks  to  the  left,  as  compared  to  Model  3  spectra  with  the 
same  density  and  temperature  parameters.  This  also  means 
the  energy  density  should  increase  near  the  left  end  with 
increasing  2.  Figure  20  portrays  the  results  of  four  flu- 
ence  spectra  of  differing  material,  with  9j=l  KeV  and  the 
same  atomic  number  density,  A/»  =  1E+27  atoms/m3. 


Fig.  20.  Comparison  of  Model  4  Spectra 

Normalized  to  Unit  Area;  «,=  i  KeV, 
Nw  =  1E+27  atoms/m3 


One  striking  feature  of  these  Model  4  specta  is  evidence 
of  a  small  peak  or  hump  at  low  energies  caused  by  the  change 
in  material  and  binding  energy,  but  primarily  binding 
energy.  (This  peak  is  not  due  to  the  same  anomaly  noted  in 
Model  1  at  the  lower  end  of  the  spectrum  because  a  smaller 
temperature  mesh  does  not  move  the  peak  closer  to  £=0  >. 

The  peak  is  more  prominent  for  higher  2  because  more  inner 
electrons  bind  per  degree  change  in  temperature,  releasing 
more  energy,  than  for  lower  2.  The  exception  in  Figure  20 
is  plutonium,  and  the  reason  it  has  no  prominent  small  peak 
is  because  the  temperature  does  not  start  out  hot  enough  to 
observe  a  rise  in  the  change  in  binding  energy  to  its  great¬ 
est  value  between  9=1  and  2  KeV  as  shown  in  Figure  21, 
although  it  does  rise  to  its  greatest  change  in  the  number 
of  binding  electrons  between  9=0,02  and  0.05  KeV ,  as  shown 
in  Figure  22. 
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Au^  /(Nv  A0  )  (keV/keV-atom ) 


Plutonium 
IE +27  atoms/m 3 


Temperature  (KeV) 

Linear-Lag  Plot  cf  Change  in  Binding  Energy 
per  Atom  with  Respect  to  Temperature 
versus  Temperature  for  Plutonium 


Anf  /(Nv  A8  )  (1 /KeV- atom) 


Plutonium 
IE +27  atome/m3 


Temperature  (KeV) 

Linear-Log  Plot  of  Change  in  Free  Electrons 
per  Atom  with  Respect  to  Temperature 
versus  Temperature  for  Plutonium 


Spectra  having  initial  temperatures  much  higher  than  1 
KeV  show  increasing  tendency  to  become  like  Model  2  spectra 


because  of  the  dominant  effect  of  radiation  energy  at  high 
temperatures  where  most  of  the  energy  is  radiated.  Spectra 
having  initial  temperatures  below  about  0.01  KeV  became  more 
like  Model  1  spectra  with  constant  material  energy  because 
all  of  the  electrons  are  bound,  eliminating  any  change  in 
binding  energy,  and  the  change  in  radiation  energy  is  too 
small  compared  to  the  change  in  material  energy. 


Again,  the  fit  is  not  as  good  as  higher  temperature  or  lower 
2  spectra  that  look  more  like  Model  2  spectra  having  only 
internal  radiation  energy.  The  fit  is  better  for  a  particu¬ 
lar  2  material  when  its  spectrum  resembles  the  plutonium 
spectrum  in  Figure  20;  that  is,  the  initial  temperature  is 
equal  to  or  lower  than  that  at  which  the  greatest  change  in 
binding  energy  occurs,  with  the  density  being  sufficiently 
low,  so  that  the  lower  end  of  the  spectrum  has  no  prominent 
small  peak  and  approaches  an  energy  density  value  close  to 
zero  . 

A  much  better  curve-fit  is  obtained  by  combining  three 
Planckians  and  a  Model  1  spectrum.  Two  Planckians  are  used 
to  match  the  curve  from  the  top  of  the  main  peak  to  the 
upper  end  of  the  spectrum,  while  a  scaled  Model  1  spectrum 
is  used  to  match  the  majority  of  the  curve  from  the  lower 
end  of  the  spectrum  to  the  top  of  the  main  peak.  The  third 
Planckian  is  scaled  to  match  the  small  peak.  The  result  is 
shown  in  Figure  24  and  compared  to  the  Model  4  copper  spec¬ 
trum  in  Figure  25. 


57 


Summary 


Model  4  presents  the  most  realistic  case  for  the  cooling 
Planckian  emitter  of  the  four  models  examined.  It  accounts 
for  losses  in  material  energy,  radiation  energy,  and  binding 
energy,  the  major  contributors  to  internal  energy,  as  the 
temperature  falls.  The  cooling  behavior  of  Model  4  may  not 
be  truly  accurate  because  of  the  somewhat  arbitrary  smooth¬ 
ing  of  the  free  electron  curves;  however,  it  gives  general 
trends  in  the  spectra  for  changes  in  temperature,  density, 
and  material  composition  parameters.  As  the  ionization 
state  falls  within  the  cooling  plasma,  the  change  in  mate¬ 
rial  and  binding  energy  may  produce  a  smaller  peak  in  energy 
density  at  lower  photon  energies  than  the  main  peak  caused 
primarily  by  the  change  in  internal  radiation  energy, 
depending  on  initial  temperature  and  density.  If  the  ini¬ 
tial  temperature  and  density  can  be  deduced  and  are  in  an 
appropriate  range,  the  material  composition  might  be 
recognized  from  its  signature  on  the  fluence  spectrum,  pro¬ 
vided  a  detector  has  the  resolution  to  d i f ferent i ate  the 
peaks . 
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VIII.  Cone lusions  and  Recommend a t ions 


The  objective  of  investigating  the  nature  of  the  time- 
integrated  spectrum  for  a  cooling  Planckian  emitter  was 
approached  by  considering  four  models  that  dealt  with  the 
composition  and  change  in  internal  energy  radiated  out 
through  the  surface  of  the  emitter.  The  change  in  material, 
radiation,  and  binding  energy  taken  separately  or  in  combi¬ 
nation,  give  insight  to  how  these  quantities  determine  the 
heat  capacity,  which  in  turn  determines  how  much  energy  is 
released  as  the  temperature  drops.  Three  parameters  were 
varied,  temperature,  density,  and  material  composition,  each 
of  which  changed  the  internal  energy  for  some  or  all  of  the 
models.  Also,  a  few  selected  spectra  were  compared  with 
Planckian  spectra  to  see  if  they  could  be  fit  with  one  or 
two  Planckian  basis  functions. 

Cone  1  us i ons 

This  study  of  t i me- i n t egr a  ted  spectra  revealed  the  fol¬ 
lowing  conclusions: 

-  The  most  significant  effect  of  material  and  binding 
energy  is  seen  at  lower  photon  energies  of  about  1  KeV 
or  less,  depending  on  initial  temperature  and  density, 
and  may  produce  a  small  peak  or  hump  in  the  spectrum 
separate  from  the  main  peak.  The  material  composition 
could  be  determined  from  this  signature  if  the  initial 
temperature  and  density  are  known. 
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-  The  fluence  spectrum  of  a  cooling  Planckian  emitter 
has  a  non-Planckian  distribution  of  energy. 

-  Two  Planckian  basis  functions  fit  fluence  spectra 
from  high  initial  temperature  emitters  best,  or  those 
that  resemble  Model  £  spectra,  whereas  the  fit  becomes 
poor  for  fluence  spectra  from  high  density  emitters,  or 
those  that  resemble  Model  1  Spectra.  In  such  case,  the 
fit  becomes  better  if  more  Planckian  basis  functions  are 
used  . 

Recommend at  ions 

Further  improvements  to  Model  4  can  be  made  by  including 
the  effect  of  expanding  volume  on  temperature  and  density, 
and  also  accounting  for  diffusion  of  energy  from  a  hotter 
interior  to  a  colder  surface  to  approximate  the  true  temper¬ 
ature  at  which  radiation  is  emitted. 
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Appendix  A:  Computer  Program  to  Calculate 
the  T i me- Integrated  Spectrum 

The  following  program  calculates  the  time-integrated 
spectrum  using  the  academic  version  of  TK  Solver  Plus  (5). 
The  heat  capacity  is  changed  in  the  subroutine  labeled  RULE 
FUNCTION:  C  for  the  different  models.  A  "C"  at  the  far 
left  edge  indicates  the  equation  is  cancelled  and  is  not 
used  in  computations,  similar  in  effect  to  a  comment  line 
denoted  with  "  . 

-  VARIABLE  SHEET  - - —  For  Acadeaic  Use  Only 

St  Input - Naae —  Output —  Unit - Coaaent - 


94 

2 

Atoaic  nuaber. 

.02 

9i 

KeV 

Initial  teaperature. 

.00001 

IohE 

KeV 

Lowest  photon  energy. 

.01 

highE 

KeV 

Highest  photon  energy. 

500 

Osteps 

Nuaber  of  teaperature  steps  (aesh). 

100 

Esteps 

Nuaber  of  energy  steps  (aesh). 

(Each  aesh  sust  be  an  even  nuaber  for 

Siapson's  rule  of  integration). 

'Planck 

func 

Planck  function  on  Function  Sheet. 

'Ionize 

func2 

Ionize  function  on  Function  Sheet. 

50 

ratio 

E/9. 

1E27 

Nv 

l/aA3 

Nuaber  density. 

area 

3.941E-8  KeV/a'3 

Area  under  non-noraalized  spectrua. 

3 . 94 10443360951 32E-8 

5.896E- 

10  areal 

KeV/a'3 

Area  1  under  non-noraalized  spectrua 
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for  split  energy  aesh. 


•ult 


RULE  SHEET 


5 . 89573572 1 9E784EE- 1 0 
Density  aultiplier. 

- . .  —  ==  For  Acadeaic  Use  Only 


S  Rule- 


Call  blank! 1 f )  "blank  erases  array  values. 

Call  blank! 'El 
Call  blank! ' s) 

Call  Cool(func;area) 

==============  PROCEDURE  FUNCTION:  Cool 

Coaaent:  Calculates  spectrua. 

Paraaeter  Variables:  lo«E,highE, Esteps, 9i,8steps 
Input  Variables:  func 

Output  Variables:  F 

S  Stateaent - 


For  Acadeaic  Use  Only 


Additional  notation  not  included  in  VARIABLE  SHEET: 
deltaE  Energy  interval. 

E  Photon  energy  value. 

’ f C i 3  (An  array).  Value  of  the  noraalized  Planck  function 

tiaes  heat  capacity  integrated  over  the  appropriate  teap- 
erature  liaits,  with  the  ith  E  as  an  energy  paraaeter. 

F  Sua  of  all  'fCil  integrated  over  the  appropriate  energy 
liaits.  (Area  under  non-noraalized  spectrua). 
k  Multiplication  factor  =  1  or  S. 

'sCil  (An  array).  Noraalized  spectrua  value  =  'fCil  /  F. 
value  Interaediate  function  value. 
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1  Check  to  see  that  energy  tesh  has  even  nuaber  of  intervals. 


IF  H0D(Esteps,2)<>0  THEN  error:=  odd_nueber_of_intervals 


1  Integration  over  energy  using  Stepson' s  rule.  Basic  procedure  is  borrowed 
‘  and  eodified  froe  reference  !5). 


del taE  :=  (highE  -  lowE)  /  Esteps 
E  :=  lowE 
k  :=  1 
value  i=  0 

FOR  i  :=  2  TO  Esteps 

E  :=  E  +  del taE 

k  ;=  3  -  k 

'fill  :=  Siepsonlfunc,  9i,  Ssteps,  E) 
value  :=  value  +  k  ♦  1 f E i ] 

NEXT  i 

■fCU  :=  Siepsonlfunc,  6i,  Ssteps,  lowE) 

1 f CEsteps  *  11  :=  Siepsonlfunc,  9i,  Ssteps,  highE) 

F  :=  (2  *  value  +  'ftll  +  'f [Esteps  ♦  ID  »  deltaE  /  3 
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Calculate  normalized  spectrue  and  fill  lists  (arrays)  for  plotting. 


E  :=  IohE 

FOR  i  :=  2  TO  Esteps  ♦  1 
E  :=  E  ♦  deltaE 
■Eli]  :=  E 
'sCil  :=  'ftil  /  F 
NEXT  i 

'ECU  :*  IohE 
'sill  :=  'fill  /  F 


Co  went: 

Parameter  Variables:  ratio, #i 
Input  Variables:  func,high#,#steps,E 

Output  Variables:  value 

S  Statement - 


PROCEDURE  FUNCTION:  Simpson  For  Academic  Use  Only 

Integration  over  temperature  using  Simpson's  rule. 


1  Description:  Standard  method  of  numerical  integration  using  second-degree 
*  polynomial  approximation  of  integrand. 


a  Additional  notation  not  included  in  VARIABLE  SHEET. 

'  apply(func,  a,  b)  Calls  a  named  function  (func)  defining  the 

'  integrand  and  supplies  the  appropriate  arguments 

■  (a,  b). 

'  delta#  Temperature  interval. 

1  high#  Upper  temperature  integration  limit. 
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r 


low©  loner  teaperature  integration  liait. 

teap  Teaperature. 

value  Intermediate  function  value  and  final  value  for  inte¬ 

gration  over  teaperature. 


‘  Check  to  see  that  teaperature  aesh  has  even  nuaber  of  intervals. 


IF  H0D(9steps,2)<>0  THEN  error:=  odd.nuaber.of.intervals 


'  Change  loner  teaperature  integration  liait  based  on  energy  E  and  ratio  of 

■  P/A  P/ratin  Biiet  ha  Inuor  frhan  Ai 


lon8  ,s  E  /  ratio 

IF  lone  >=  9i  THEN  error  :=  ratio. is.too.saall 


*  Integration  over  teaperature  using  Simpson ' s  rule.  Basic  procedure  is  bor 

•  roned  and  modified  froa  reference  !5). 


delta0  :=  (high9  -  lon8)  /  9steps 
teap  :=  lon0 
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k  :=  1 


value  :=  0 

FOR  i  :=  2  TO  Ssteps 
teap  :=  teap  +  del  tad 
k  :=  3  -  k 

value  :=  value  +  l  »  applylfunc,  E,  teap,  deltas,  lo*S,  highS) 

NEXT  i 

value  :=  (2  *  value  ♦  applylfunc,  E,  load,  deltas,  lowS,  highS)  +  applylfunc, 

'  E,  highS,  deltas,  lowS,  highS) >  *  deltas  /  3 

.  . ■  RULE  FUNCTION:  Planck - For  Acadeaic  Use  Only 

Coaaent:  Noraalized  Planck  function. 

Paraaeter  Variables: 

flrguaent  Variables:  E, S, deltas, lowS, highS 
Result  Variables:  p 

S  Rule - 

'  Additional  notation  not  included  in  VARIABLE  SHEET: 

1  p  Noraalized  Planckian  tiaes  heat  capacity. 

'  planck  Noraalized  Planckian  value. 


planck  8  15  /  (piOPA  *  E“3  /  <S‘A  ♦  (EXPIE/S)  -  1)> 


p  8  planck  a  CIS,  deltas,  IomS,  highS) 


RULE  FUNCTION:  C 


Coaaent: 

Paraaeter  Variables: 
Arguaent  Variables: 
Result  Variables: 


Heat  capacity. 
Nv,Si,Z 

S, deltas, IohS, highS 
heatcap 


For  Acadeaic  Use  Only 
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Rule 


Additional  notation  not  included  in  VARIABLE  SHEET  and  previous  functions. 
a£  Radiation  constant, 

heatcap  Heat  capacity, 

ionl  Binding  energy  at  teapl. 

ionS  Binding  energy  at  teapE. 

Ionize  Function  that  computes  binding  energy  per  atoi. 

•ateriall  Material  energy  density  at  teepl. 

aaterialS  Material  energy  density  at  teapE. 

radl  Radiation  energy  density  at  teapl. 

rad£  Radiation  energy  density  at  teapE. 

starA  Interaediate  Zstar  function  value. 

starB  Interaediate  Zstar  function  value. 

starl  Zstar  function  value  for  teapl. 

starS  Zstar  function  value  for  teapE. 

teap  Interaediate  teaoerature  value, 

teapl  Lower  teaperature  of  teaperature  interval. 

teapE  Upper  teaperature  of  teaperature  interval. 

teap3  Interaediate  teaperature  value. 

9  Teaperature  value. 

Zstar  Function  that  coaputes  the  nuaber  of  free  electrons 
per  atoa  at  a  given  teaperature. 


For  Model  1. 
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C  heatcap  *  1 


1  For  Model  2. 

heatcap  s  <T3 


*  For  Model  3. 


C  a2  =  8.5634E+2B  'Radiation  constant  expressed  in  l\(KeV-t)‘3. 
C  heatcap  *  Nv  *  (3  /  2)  ♦  4  *  a2  »  0“3 

'  For  Model  4. 


C  a2  *  8.5634E+28  'Radiation  constant  expressed  in  l\(KeV-e)'3. 


'  Detereine  lower  and  upper  temperatures  of  interval. 


C  teap  =  8  -  del  tad  /  2 


69 


C  IF  9  =  load  THEM  teapl  =  load  ELSE  tetpl  -  teap 
C  teap 3  =  d  +  del  tad  /  2 

C  IF  d  »  highd  THEN  teap2  =  highd  ELSE  teapH  =  teap3 


1  Deteraine  nuaber  of  free  electrons  per  atoa  at  both  teaperatures. 


C  star A  =  A8S(Zstar( teapl)) 

C  IF  star A  >  Z  THEN  starl  =  Z  ELSE  starl  =  star A 
C  starB  =  ABS(Zstar(teap2H 
C  IF  starB  >  Z  THEN  star2  =  Z  ELSE  star2  =  starB 


1  Coapute  aaterial  energy  density  at  both  teaperatures. 


C  aatterl  =  Nv  *  (1  +  starl)  »  3  /  2  *  teapl 
C  «atter2  *  Nv  »  (1  ♦  star2)  t  3  /  2  «  teap2 


'  Coapute  radiation  energy  density  at  both  teaperatures. 


C  radl  -  a2  *  teapl" 4 
C  rad2  1  a2  *  teap2"4 
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Compute  binding  energy  per  atoa  at  both  temperatures. 


C  ioni  =  Nv  *  ABSdonize(starli) 
C  ion2  =  Nv  *  ABS(Ionize(star2) ) 


*  Coapute  heat  capacity  for  temperature  interval. 

■  ssssssssssssssssssssssssssssssssssssssssssrssss 


C  heatcap  =  (aatter2  *  rad2  ♦  ion2  -  (matter!  *  rad!  +  ioni ) )  /  del ta6 
. . .  „a.m ,  -rrawBmr — t  LIST  FUNCTION:  2star  _  For  Academic  Use  Only 


Comment: 

Interpolation  function  for  free  electrons  vs  temperature. 

Domain  List: 

kT7a 

Napping: 

Linear 

Range  List: 

nf7a 

Element—  Domain — 

-  Range - 

1  0 

0 

2  .00245345426244727  .01 

3  .0040090803003785  .1 

4  .00766489707885056  .5 

5  .0113754622168454  1.5 

6  .0133180652203634  2 

7  .0168206934563477  3 

8  .0240665353092386  5 
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' 

9 

.0475779982016186 

10 

• 

10 

.0857295777914731 

15 

11 

.113199665581666 

17 

la 

.177949883904736 

19.5 

13 

.278790992266599 

23 

14 

.348411231622957 

25 

15 

.407371098297251 

26 

16 

.465070257418233 

26.5 

17 

1.36937648331806 

28.5 

18 

1.70020947649073 

28.9 

19 

2.12480361177458 

28.98 

20 

2.36418313794776 

28.99 

21 

3.62590540771914 

29 

• 

22 

10 

29 

.  - = 

—  . LIST  FUNCTION:  Ionize  s=s==  For  Acadeaic  Use  Only  1 

Coaaent:  Interpolation  function  for  binding  energy  vs  nf. 

Doaain  List:  nfCu 

Happing:  Linear 

Range  List:  bindCu 

Element— 

Doaain - 

-  Range - 

1 

0 

0 

2 

1 

.01055774 

3 

2 

.10225568 

4 

3 

.20641458 

5 

4 

.32435908 

• 

6 

5 

.45741488 

72 

7 

6 

.60690848 

8 

7 

.77416718 

9 

8 

.96051848 

10 

9 

1.16729018 

11 

10 

1.39581038 

12 

11 

1.64740738 

13 

12 

1.92340938 

14 

13 

2.22514478 

15 

14 

2.55394218 

16 

15 

2.91112998 

17 

16 

3.29803688 

18 

17 

3.71599148 

19 

18 

4.16632238 

20 

19 

4.65035818 

21 

20 

6.26254918 

22 

21 

7.96334218 

23 

22 

9.75521318 

24 

23 

11.64063918 

25 

24 

13.62209718 

2  4 

25 

15.70206518 

27 

26 

17.88301918 

28 

27 

20.16743618 

29 

28 

31.01259618 

30 

29 

42.09627618 
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Appendix  B: 


with  T  wo  P 1 anc  k i an  Bas i s  Func  t ions 

This  program  computes  the  best  least  squares  approxima¬ 
tion  for  a  spectrum  from  a  given  range  of  temperatures  and 
two  Planckian  basis  functions.  Its  output  gives  the  two 
Planckian  temperatures  and  the  coefficients  necessary  to 
calculate  the  curve,  plus  a  measure  of  how  good  the  fit  is 
using  the  sum  of  the  square  of  the  errors  or  difference 
between  original  and  calculated  Y-axis  values.  The  program 
is  written  in  FORTRAN  77  language  and  calls  a  mathematical 
library  routine  DFNLSQ  written  by  IMSL  (6).  Typical  output 
follows  the  program. 

*  Maj  Drew  Fisher 

*  14  Oct  89 
PROGRAM  PLANCK 

*  ============ ================3=====-============ ================ 

*  This  program  computes  a  least  squares  approx imat ion  with  user- 

*  supplied  basis  functions.  The  approximation  takes  the  form 

* 

*  Y  =  COEFF(l)  +  COEFF (2 ) *F ( 1 , X )  +  COEFF ( 3 ) «F ( B , X ) 

*  Definitions: 


* 

COEFF 

Coefficient  array. 

# 

C0EFF2 

Storage  array  for  coefficients. 

* 

DELTAT 

Temperature  increment. 
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4 

DFNLSQ 

— 

Library  routine  that  computes  least 

4 

squares  approximation  with  user-supplied 

* 

basis  functions. 

4 

ENERGY 

- 

Array  of  X-axis  values. 

* 

F 

- 

Subroutine  for  computing  values  from 

# 

basis  functions. 

# 

FRACT 

- 

Array  of  original  Y-axis  values. 

4 

INDEX  1 

- 

Stores  index  counter  I. 

4 

INDEX2 

- 

Stores  index  counter  J. 

4 

INTCEP 

- 

Intercept . 

4 

IWT 

- 

Weight  option  (not  used). 

4 

JJ 

- 

A  counter. 

4 

N 

- 

Counter  used  for  number  of  data  points. 

* 

NBAS IS 

- 

Number  of  basis  functions. 

* 

NFL  IP 

- 

Flip  flop  counter. 

4 

NSIZE1 

- 

Array  size. 

4 

NS  I ZEE 

- 

Array  size. 

4 

NUM1 

- 

Loop  index. 

4 

NUM2 

- 

Loop  index. 

4 

RESET 

- 

Original  TEMPI. 

* 

SSE 

- 

Sum  of  squares  of  the  errors  -  an  indicator 

4 

of  fit. 

* 

STORE 

- 

Storage  array  for  SSE. 

4 

TEMPI 

- 

Stores  temperature  for  first  basis  function. 

4 

TEMPE 

- 

Stores  temperature  for  second  basis  function 

4 

THETA  1 

- 

Temperature  for  first  basis  function. 
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THETA2 


Temperature  for  second  basis  function 


* 

*  WEIGHT  -  Weights  given  to  basis  functions  (an  array  of 

*  higher  dimension  if  weights  other  than  1  are 

*  assigned ) . 


INTEGER  I,  INDEX  1 ,  INDEX2,  INTCEP,  IWT,  J,  JJ,  M,  N,  NBASIS 
INTEGER  NFL  IP,  NSIZE1 ,  NSIZE2 ,  NUMI ,  NUM2 

PARAMETER  (IWT  =  0,  NBASIS  =  2,  NSIZE1  =  10,  NSIZE2  =  200) 
DOUBLE  PRECISION  C0EFF(NSIZE1 ) ,  CQEFF2(NSIZEl ) ,  DELTAT 
DOUBLE  PRECISION  ENERGY ( NSIZE2 ) ,  FRACT ( NSIZE2 ) ,  RESET,  SSE 
DOUBLE  PRECISION  ST0RE(2),  TEMPI,  TEMP2 ,  THETA 1,  THETA2 ,  WEIGHT 
PARAMETER  (WEIGHT  =  1) 

COMMON  THETA 1,  THETA2 
EXTERNAL  DFNLSQ,  F 


*  Open  file  containing  X-axis  and  Y-axis  spectrum  values. 


OPEN  (UNIT  =  20,  FILE  =  ‘cufit‘,  STATUS  =  ‘OLD’) 

N  =  1 

10  READ  (UNIT  =  20,  *,  END  =  20)  ENERGY ( N ) ,  FRACT ( N ) 
PRINT*,  N,  ENERGY <N) ,  FRACT ( N ) 

N  =  N  +  1 
GOTO  10 
20  CONTINUE 
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N  =  N  -  1 


PRINT* 

PRINT*,  'Number  of  points  =  ',  N 
CLOSE  (UNIT  =  SO) 


User  input. 


PRINT*,  ‘Include  intercept?  (1  =  Yes,  0  =  No)' 

PRINT* 

READ*,  INTCEP 
PRINT* 

PRINT*,  'For  investigation,  start  THETA1  and  THETAE  at' 
PRINT*,  'what  temperatures  (KeV)?  ( THETA 1 ,  THETAE) ' 

PRINT* 

READ*,  THETA 1,  THETAS 
PRINT* 

PRINT*,  'Increment  temperatures  by  what  DELTAT?' 

PRINT* 

READ*,  DELTAT 
PRINT* 

PRINT*,  'Loop  information.  How  many  iterations  for  each' 
PRINT*,  'loop?  (inside  C THETA 1 1 ,  outside  C THETAS] ) ' 

PRINT* 

READ*,  NUMS,  NUM1 
PRINT* 
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Initialize  values 


* 

* 

* 


NFL IP  = 

1 

RESET  = 

THETA 1 

STORE! 1 ) 

=  l.OD+OO 

STORE (2) 

=  STORE! 1) 

TEMPI  = 

THETA 1 

TEMP2  = 

THETAS 

INDEX  1  = 

1 

INDEX2  = 

1 

JJ  =  0 

*  Start  search  for  temperatures  that  give  lowest  SSE. 


t 


DO  00  I  =  1,  NUM1 
PRINT*,  'I  =  '  ,  I 
THETA 1  =  RESET 
DO  70  J  =  1,  NUM2 

IF  ( THETA 1  .NE.  THETA2 )  THEN 
JJ  =  JJ  +  1 

CALL  DFNLSQ  (F,  INTCEP,  NBASIS, 
WEIGHT,  COEFF,  SSE) 


N, 


ENERGY,  FRACT, 


I WT , 


70 


Store  coefficients  from  first  loop  pass  in 
temporary  storage. 


IF  (JJ  .EQ.  1)  THEN 

DO  40  M  =  1,  NBASIS  +  1 
C0EFF2(M)  =  COEFFCM) 
CONTINUE 
ENDIF 

NFL IP  =  3  -  NFL IP 
STORE! NFL  IP)  =  SSE 


Compare  SSE  with  the  old  value  of  SSE  and  overwrite 
temporary  storage  locations  with  the  smallest  value 
and  its  associated  values.  Three  cases  are  checked 
and  a  fourth  one  implies  no  change. 


IF  (NFL IP  .EQ.  2  .AND.  ST0RE(2)  .LT.  STORE ( 1 ) )  THEN 
STORE < 1 )  =  STORE (2) 

TEMPI  =  THETA  1 
TEMP2  =  THETAS 
INDEX  1  =  I 
INDEX2  =  J 


DO  50  M  =  1,  NBASIS  ♦  1 


COEFFS(M)  =  COEFF(M) 


50  CONTINUE 

ELSE IF  (NFL IP  .EQ.  1  .AND.  STORE(l)  .LT.  STORE ( 2 )  )  THEN 
TEMPI  =  THETA 1 
TEMP2  =  THETA2 
INDEX  1  =  I 
INDEX2  =  J 

DO  60  M  =  1,  NBAS IS  +  I 
C0EFF2(M)  =  COEFF(M) 

60  CONTINUE 

ELSE IF  (NFL IP  .EQ.  1  .AND.  ST0RE(2>  .LT.  STORE ( 1 ) )  THEN 
STORE ( 1 )  =  ST0RE(2) 

ELSE 
END  IF 

STORE (NFL  IP)  =  STORE(l) 

END  IF 

*  ========================================================= 

*  Increment  temperature  for  inner  loop,  and  after  last  pass 

*  print  smallest  SSE  encountered  in  inner  loop,  then 

*  increment  temperature  for  outer  loop. 

*  ========================================================= 

THETA1  =  THETA 1  +  DELTAT 
70  CONTINUE 

PRINT#,  STORE ( 1 ) 

THETA2  =  THETA2  +  DELTAT 
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80  CONTINUE 


*  ============= 

*  Print  output. 

*  ============= 


PRINT* 

PRINT*, 

'I  =' ,  INDEX  1 ,  ' J 

=‘,  INDEX2 

PRINT*, 

' THETA 1  =',  TEMPI, 

' THETA2  =' , 

TEMPS 

PRINT*, 

‘Sum  of  the  square 

of  errors  = 

1  ,  STORE ( 1 ) 

PRINT*, 

'Coefficients  are: 

• 

DO  90  M 

=  1,  NBAS IS  +  1 

PRINT#,  C0EFF2<  M) 

90  CONTINUE 
END 

*  =============================== 

*  Subroutine  for  basis  functions. 

* 


* 

* 

* 

# 

* 

* 

# 


Definitions: 

F  -  Subroutine  for  computing  values  from  user 

supplied  basis  functions. 

K  -  Basis  function  number. 

X  -  X-axis  value. 

THETA1  -  Temperature  for  first  basis  function. 
THETAS  -  Temperature  for  second  basis  function. 
DOUBLE  PRECISION  FUNCTION  F(K,  X) 
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INTEGER  K 


DOUBLE  PRECISION  THETA 1,  THETA2,  X 
COMMON  THETA 1,  THETAS 
IF  (K  .EQ.  1)  THEN 

F  =  X**3.0D+00  /  ( THETA l**4.0D+00  *  <DEXP( X/THETA1 ) 

$  -  l.OD+OO)) 

ELSE 

F  a  X**3.0D+00  /  < THETA2**4 . OD+OO  *  (DEXP( X/THETA2) 

*  -  1 .OD+OO) ) 

END  IF 
RETURN 
END 

The  following  output  is  typical  of  the  above  program  and 
was  used  to  generate  the  two-Pl anck i an  fit  in  Figure  10  of 
Chapter  V. 

*/.  f77  planck.f  -limsl 
p 1 anck . f : 

MAIN  p 1 anck : 
f : 

*/.  a.  out 

Number  of  points  =  197 

Include  intercept?  <1  =  Yes,  0  =  No) 

0 

For  investigation,  start  THETA1  and  THETA2  at 
what  temperatures  (KeV)?  (THETAl,  THETA2) 
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.409,  .840 


Increment  temperatures  by  what  DELTAT? 

.001 

Loop  information.  How  many  iterations  for  each 
loop?  (inside  CTHETA1D,  outside  CTHETA23) 

3 ,  6 
I  =  1 

3 . 09545336704 1 7d-04 
I  =  2 

3 . 00249629 1 6826d-04 
I  =  3 

2 . 9380294559385d-04 
I  =  4 

2. 901 9676931 35 ld-04 
I  =  5 

2 . 8935 1 58596325d— 04 
I  =  6 

2 . 8935 1 58596325d-04 
I  =  5  J  =  2 

THETA  1  =  0.41000000000000  THETA2  =  0.84400000000000 

Sum  of  the  square  of  errors  =  2 . 8935 1 58596325d-04 

Coefficients  are: 

1 . 8698923739682d-02 
0.  13468654781493 
0. 

The  two-P 1 anck i an  fit  is  calculated  from: 
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0 


a0 

CEl 


a2 

6i 

02 

Pi 

Pz 

Sm(f) 


1 . 8698923739682 F  -  02 

0.13468654781493 

0.410 

0.844 

FVtexpCF/e^-l] 
F3/[exp(F/02)-  1] 
aQ  +  a,  F,  +  a2  *  P2 
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Appendix  C 


Energy  and  Temperature 

The  following  computer  program  calculates  the  binding 
energy  to  remove  the  next  bound  electron  from  an  atom  and  is 
written  in  Quick  Basic  (7). 

'  Maj  Drew  Fisher 

*  5  Nov  89 

*  2STAR.BAS 

'  Ver .1.0  5  Nov  89 

'  This  program  computes  Zj*  and  Ej  for  use  in  Zstar.tk  to 
find  free 

'  electrons  vs  temperature. 

OPTION  BASE  1 
DEFINT  I  — N 
SIZE  =  10 

REDIM  xN(SIZE),  SIGMACSIZE,  SIZE) 


1  Slater  screening  constants  from  (3:169). 
'  Column  1-5. 


DATA 

0.6250, 

0.9383, 

0.9811, 

0.987  , 

0.994 

DATA 

0.2346, 

0.6895, 

0.8932, 

0.94  , 

0.97 

DATA 

0. 1090, 

0.3970, 

0.7018, 

0.85  , 

0.92 

DA  r  A 

0.0617, 

0.2350, 

0.4781  , 

0.705  , 

0.83 

DATA 

0.0398, 

0.  1552, 

0.3312, 

O .531  , 

0.72 

DATA 

0.0277, 

0. 1093, 

0.2388, 

0.400  , 

0.854 
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DATA 

0.0804, 

0.0808, 

0. 1788, 

0.3108, 

0.459 

DATA 

0.0156, 

0.0685, 

0. 1378, 

0.8485, 

0.371 

DATA 

0.0183, 

0.0494, 

0. 1  106, 

0. 1936, 

0.899 

DATA 

0.0100, 

0.0400, 

0.0900, 

0.1584, 

0.845 

'  Column  6  - 

10. 

DATA 

0.997  , 

0.999  , 

1.000  , 

1.000  , 

1  .000 

DATA 

0.984  , 

0.990  , 

0.993  , 

0.995  , 

1 .00 

DATA 

0.955  , 

0.97  , 

0.98  , 

0.99  , 

1 .00 

DATA 

0.90  , 

0.95  , 

0.97  , 

0.98  , 

0.99 

DATA 

0.83  , 

0.90  , 

0.95  , 

0.97  , 

0.98 

DATA 

0.735  , 

0.83  , 

0.90  , 

0.95  , 

0.97 

DATA 

0.610  , 

0.745  , 

0.83  , 

0.90  , 

0.95 

DATA 

0.506  , 

0.635  , 

0.750  , 

0.83  , 

0.90 

DATA 

0.431  , 

0.544  , 

0.656  , 

0.760  , 

0.83 

DATA 

0.353  , 

0.466  , 

0.576  , 

0.67  , 

0.765 

'  Read  data  into  array  SIGMAil,  J). 
CLS 

FOR  I  =  1  TO  10 

FOR  J  =  1  TO  5 

READ  S I GMA ( I  ,  J) 

NEXT  J 
NEXT  I 

FOR  I  =  1  TO  10 

FOR  J  =  6  TO  10 
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READ  S I GMA < I ,  J) 


NEXT  J 
NEXT  I 


'  User  input. 

INPUT  "Atomic  number,  Z" ;  Z 
PRINT 

INPUT  "Level,  J";  J 
PRINT 

INPUT  "Number  of  free  electrons  per  atom";  xNF 
PRINT 

INPUT  "Number  density  Cl/m~33";  xNUMDEN 
PRINT 


'  Compute  effective  atomic  number  Zj*. 

sum  =  0 

FOR  1=1  TO  J 

xN ( I )  =2  *  I  ~  2 

TEMP  =  SUM  1 

SUM 1  =  SUM  1  +  xN( I ) 

IF  Z  -  xNF  <  SUM  1  THEN  xN(I)  =  Z  -  xNF  -  TEMP 
NEXT  I 
SUMS  =  0 
IF  J  >  1  THEN 

FOR  I  =  1  TO  J  -  1 
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SUMS  =  SUMS  +  SIGMAd,  J>  *  xN(I) 

NEXT  I 
END  IF 

2 j STAR  =  Z  -  SUM2  -  xN(J)  *  < 1  -  1  /  < 2  *  J  ~  2 ) )  *  SIGMA(J, 
J) 


’  Compute  ionization  potential,  Ej,  to  remove  electron. 

aO  =  5.292E-11  'Bohr  radius  Cm3. 

eCHARGE  =  1.602E-19  'Electron  charge  EC3  . 

EPSILO  =  8.854E-12  'Permittivity  of  vacuum  CF/m] 

PI  =  3.14159 

RYDBERG  =  13.61  'Rydberg  constant  CeV] . 

rj2  =  aO  ~  2  *  J  ~  4  /  2 j STAR  ~  2  *  (7/4+5/  ( 4  *  J  ~ 
2)  ) 

R  =  <3/  <4  *  PI  *  xNUMDEN ) )  ~  (1/3) 

ENER  =  eCHARGE  /  (8  *  PI  *  EPSILO  *  R)  'CeV3 

Ej  =  -RYDBERG  *  ( 2 jSTAR  /  J )  ~  2  +  xNF  *  ENER  *  <18  /  5  - 

r  j  2  /  R  ~  2) 


'  Print  output. 

PRINT  "Atomic  number  =  ";  2 

PRINT  "Number  density  Cl/m~33  =  xNUMDEN 
PRINT  "Level,  J  =  J 

PRINT  "Electrons  in  J  level  <nj>  =  xN(J) 
PRINT  "Free  electrons  per  atom  (2*)  =  xNF 
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PRINT  "Effective  atomic  number  (Zj*)  =  ZjSTAR 
PRINT  "Ionization  potential  <Ej>  CeV]  =  ";  Ej 
PRINT 
END 


The  next  program  computes  the  temperature  corresponding 
to  the  binding  energy  calculated  in  the  previous  program. 


Pt  Innut _ 

_  Naaa _ 

VARIABLE  SHEET  = 

Output —  Unit 

■  ■■  - . . For  Acadeaic  Use  Only 

5.68SE-15  »e 

see  >> 

Electron  aass  in  KeV(s/«)"2. 

4.136E-18  h 

KeV-s 

Planck  constant. 

29 

z 

Atoaic  nuaber. 

-10860.9 

Ej 

eV 

Ionization  potential. 

L  27.1 

nf 

1/atoa 

Free  electrons  per  atoa. 

1EH7 

Nv 

l/a“3 

Atoaic  nuaber  density. 

1 

j 

Principal  quantua  state. 

1.9 

nj 

Nuaber  of  electrons  in  state  i. 

L 

9 

.92992099  KeV 

Teaperature. 

P 

.00038227 

Parameter. 

Rill  F  SWFFT 

;  Pulp 

P  =  h“2  / 

(2  *  pi!) 

*  ae)  *  (nf  »  Nv 

/  2  •  (2  *  j"2  /  nj  -  1))A(2  /  3) 

Q  =  -Ej  /  ln(ie  /  Pr(1.5)) 
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